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1 Introduction 

A quantum mechanical system of N particles in d dimensions can be described by a complex valued 
wave function ijj^ € L 2 (M. dN , dx\ . . . dx^). The variables xi, . . . , xn £ K d represent the position of 
the N particles. Physically, the absolute value squared of iPn(xi,X2, ■ ■ ■ , xjv) is interpreted as the 
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probability density for finding particle one at x\, particle two at x%, and so on. Because of this 
probabilistic interpretation, we will always consider wave functions ipw with L 2 -norm equal to one. 

In Nature there exist two different types of particles; bosons and fermions. Bosonic systems are 
described by wave functions which are symmetric with respect to permutations, in the sense that 

1p N (x lTl ,X lT 2,---,X lTN ) = Tp N (xi,...,X N ) (1.1) 

for every permutation ir € SV- Fermionic systems, on the other hand, are described by antisymmetric 
wave functions satisfying 

V>7v0r7ri,£7r2,---,a;7r7v) = a^ N (xi, . . . ,x N ) for all IT G S N , 

where ov is the sign of the permutation ir; a n = +1 if ir is even (in the sense that it can be written 
as the composition of an even number of transpositions) and a n = — 1 if it is odd. In these notes 
we are only going to consider bosonic systems; the wave function ip^ will always be taken from the 
Hilbert space L^(R ), the subspace of L 2 (M ) consisting of all functions satisfying (11. ip . 

Observables of the iV-particle system are self adjoint operators over Lg(R ). The expected value 
of an observable A in a state described by the wave function ip^ is given by the inner product 

(ip N ,AtpN) = / dxi ...dx N ip N (xi,... ,x N ) (Aip N )(x 1 , . . . ,x N ). 

The multiplication operator Xj is the observable measuring the position of the j'-th particle. The 
differential operator pj = — zVj is the observable measuring the momentum of the j-th particle (j)j 
is called the momentum operator of the j-th particle). 

The time evolution of an iV-particle wave function ipx E L 2 (M ) is governed by the Schrodinger 
equation 

id t tpN,t = H N ip N>t . (1.2) 

Here, and in the rest of these notes, the subscript t indicates the time dependence of the wave 
function; all time-derivatives will be explicitly written as dt- On the right hand side of (|1.2p . H n is a 
self-adjoint operator acting on the Hilbert space L 2 (1R ), usually known as the Hamilton operator 
(or Hamiltonian) of the system. We will consider only time-independent Hamilton operators with 
two body interactions, which have the form 

N N 

H N = Y1 (" A ^ + ^xt^j)) + A Y, V(Xi ~ Xj) . 

j=l i<j 

The first part of the Hamiltonian is a sum of one-body operators (operators acting on one particle 
only); the sum of the Laplacians is the kinetic part of the Hamiltonian. The function V^xt describes 
an external potential which acts in the same way on all iV particles. For example, V ex t may describe a 
confining potential which traps the particles in a certain region. The second part of the Hamiltonian, 
given by a sum over all pairs of particles, describes the interactions among the particles (A £ K is 
a coupling constant). The Hamilton operator is the observable associated with the energy of the 
system. In other words, the expectation 



(i/) N ,H N tp N ) = / cLci . . . daw ip N (xi, . . . ,x N )(H N ijj N )(xi, . . . ,x N ) 
gives the energy of the system in the state described by the wave function tp^. 



Note that the Schrodinger equation (|1.2p is linear and, since Hjy is a self-adjoint operator, it 
preserves the L 2 -norm of the wave function (moreover, since Hn is invariant with respect to permu- 
tations, it also preserves the permutation symmetry of the wave function). In fact, the solution to 
(jl.2p . with initial condition f/ ; A r ,t=o = "0iV> can be written by means of the unitary group generated 
by H]\f as 

$ Nt = e~ iHNt ilj N for allieM. (1.3) 



The global well-posedness of (JL2J) is not an issue here. The study of (|1.2|> is focused, therefore, on 
other questions concerning the qualitative and quantitative behavior of the solution ipN,t- Despite 
the linearity of the equation, these questions are usually quite hard to answer, because in physically 
interesting situation the number of particles N is very large; it varies between N ~ 10 3 for very 
dilute Bose-Einstein samples, up to values of the order N ~ 10 30 in stars. For such huge values of N, 
it is of course impossible to compute the solution (|1.3|) explicitly; numerical methods are completely 
useless as well (unless the interaction among the particles is switched off). 

Fortunately, also from the point of view of physics, it is not so important to know the precise solu- 
tion to (|1.2|) : it is much more important, for physicists performing experiments, to have information 
about the macroscopic properties of the system, which describe the typical behavior of the particles, 
and result from averaging over a large number of particles. Restricting the attention to macroscopic 
quantities simplifies the study of the solution ipN,t, but it still does not make it accessible to math- 
ematical analysis. To further simplify matters, we are going to let the number of particles iV tend 
to infinity. The macroscopic properties of the system, computed in the limiting regime N — > oo, are 
then expected to be a good approximation for the macroscopic properties observed in experiments, 
where the number of particles iV is very large, but finite (in some cases it is possible to obtain explicit 
bounds on the difference between the limiting behavior as N — ► oo and the behavior for large but 
finite N; see Section l4~2l) . 



To consider the limit of large N, we are going to make use of the marginal or reduced density 
matrices associated with an N particle wave function ipN £ L^(R ). First of all, we define the 
density matrix 7^ = iV'ivKV'ivl associated with ip^ as the orthogonal projection onto ipN', we use 
here and henceforth the notation \ip){ip\ to indicate the orthogonal projection onto ip (Dirac bracket 
notation). Note that, expressed through the density matrix 77V, the expectation (ij;n,Aijjn) of the 
observable A can be written as TV A^n. The kernel of the operator 7jv is then given by 

7jv(x;x') = tp(x)tp(x') , 

where we introduced the notation x = {x\, . . . , xjy),x' = (x[, . . . ,x' N ) S ~R dN . Note that the L 2 - 
normalization of tpN implies that Tr7jy = 1. For k = 1, . . . ,N, we define the ^-particle marginal 

(k) 

density 7^ associated with ■i/'iv as the partial trace of jn over the degrees of freedom of the last 
(N — k) particles: 

7^- ) = Tr fc+ i ife+2 ,...,Ar \iI>n)(i/>n\ 

(k) 

where Trjfc-fi jv denotes the partial trace over the particle k + 1, k + 2, . . . , N. In other words, ^y N 
is defined as the non-negative trace class operator on L 2 (\R ) with kernel given by 

7y(x fe ;x' fc ) = / dxjv_ fe 7Af( x fc> x iV-fc; x fc) x A r -fc) 

(1.4) 

dxjy^fc f/>jv( x fc> *N-k) i>N{*k> x N-k) ■ 



The last equation can be considered as the definition of partial trace. Here we used the notation 
x fc = (xi,...,Xk),x! k = (x[, . . . , x' k ) G R dk and *N-k = (xk+l, ■■■ i x n) G M d(Ar_fc) . By definition, 



(k) 

Trjjy = 1 for all N and for all k = 1, . . . ,N (note that, in the physics literature, one normally 
uses a different normalization of the reduced density matrices). For fixed k < N, the /c-particle 

density matrix does not contain the full information about the state described by ipjy. Knowledge 

(k) 
of the /c-particle marginal 7^ , however, is sufficient to compute the expectation of every /c-particle 

observable in the state described by the wave function ipj^. In fact, if A^ k > denotes an arbitrary 

bounded operator on L 2 (M ), and if A^ ' ® l( N ~ k > denotes the operator on L 2 (M ) which acts as 

A^ k ' on the first k particles, and as the identity on the last (N — k) particles, we have 

(^V, (a<*> Si 1<""*)) tf, N ) = Tr (a« ® 1<"-*)) 1N = Tr A« 7 «. (LB) 

(k) 

Thus, 7^ is sufficient to compute the expectation of arbitrary observables which depend non-trivially 
on at most k particles (because of the permutation symmetry, it is not important on which particles 
it acts, just that it acts at most on k particles). 

Marginal densities play an important role in the analysis of the N — > 00 limit because, in contrast 

(k) 

to the wave function ip^ and to the density matrix 77V, the /c-particle marginal 7jy can have, for 
every fixed k £ N, a well-defined limit as N —>■ 00 (because, if we fix k £ N, {7^ } defines a sequence 
of operators all acting on the same space L 2 (R )). 

In these notes we are going to study macroscopic properties of the dynamics of bosonic iV-particle 

(k) 

systems, in the limit N — ► 00. We are interested in the time-evolution of marginal densities "y N t 
associated with the solution i/)N,t to the Schrodinger equation (|1.2j) . for fixed k, and as N — ► 00. 

(k) 

Unfortunately, it is not so simple to describe the time-dependence of jj^ t in the limit of large TV; in 
fact it is in general impossible to obtain closed equations for the evolution of the limiting /c-particle 

(k) (k) 

density 7^ t = limjv_ +00 ^ N t (in general it is not even clear that this limit exists). Nevertheless, there 

(k) 

are some physically interesting situations for which it is indeed possible to prove the existence of 7^ t 
and to derive closed equations to describe its dynamics. In Section [21 we are going to study the time 
evolution of factorized initial wave functions of the form ^at(x) = n,-_i f( x j) m the so-called mean 
field limit. We will show that in this case, for every fixed k £ N, the /c-particle marginal associated 
with the solution to the Schrodinger equation ipN,t converges, as N — > 00, to the limiting /c-particle 

density 7^ = \(ft)(<pt\ , where tpt is the solution of a certain one-particle nonlinear Schrodinger 
equation, known as the Hartree equation. In Section [3l we are going to study the time-evolution 
of Bose-Einstein condensates, in the so-called Gross-Pitaevskii limit. As we will see, although it 
describes a very different physical situation, the Gross-Pitaevskii scaling can be formally interpreted 
as a mean- field limit, with a very singular interaction potentials. Also in this case we will prove 
that the time evolution of the marginal densities can be described through a one-particle nonlinear 
Schrodinger equation (known as the time-dependent Gross-Pitaevskii equation); however, because 
of the singularity of the interaction potential, the analysis in this case is going to be much more 
involved. In Section HI we will come back to the study of mean field models, and we will discuss how 
to prove quantitative estimates on the rate of convergence towards the Hartree dynamics. 

Notation. Throughout these notes, we will make use of the notation x = (x\, ■ ■ ■ ,xn),x! = 
(xi,. . .x' N ) £ R dN , and for k = 1, . . . ,N, x ft = (xi,. . . ,x k ),x' k = (x[, . . . ,x' k ) £ R dk , and x N _ k = 
(xfc+i, . . . , xn) £ ]]£"(" -fc ) (starting from Section El we will fix d = 3). We will also use the shorthand 
notation Vj = V x . . and Aj = A x . . 



2 Derivation of the Hartree Equation in the Mean Field Limit 

2.1 Mean Field Systems 

A mean-field system is described by an iV-body Hamilton operator of the form 



N , N 



H N = J2 (-Aj + V cxt ( Xj )) + ^fl V & ~ x ^ ( 2 - X ) 



N 

j = l Kj 

acting on the Hilbert space L 2 (R ), the subspace of L 2 (K ) consisting of permutation symmetric 
functions. In these notes we will only discuss bosonic systems, which are described by symmetric wave 
functions. Note, however, that the mean field limit for fermionic system has also been considered in 
the literature; see, for example, [26j ED U\ • In (|2.1|) and henceforth, we use the notation Aj = A Xj 
(similarly, we will use the notation Vj = V x ). The mean-field character of the Hamiltonian is 
expressed by the factor 1/N in front of the interaction; this factor guarantees that the kinetic and 
potential energies are typically both of order N. 

We are going to be interested in the solution tpNj = e~ Nt t/)N of the Schrodinger equation (II. 2j) 
with Hamiltonian Hn given by (|2.ip and with factorized initial data 

N 

^v(x) = JJ (p( Xj ) , (2.2) 

i=i 

for some ip £ L 2 (R d ). The physical motivation for studying the evolution of factorized wave functions 
is that states close to the ground state of Hn (the eigenvector associated with the lowest eigenvalue) , 
which are the most accessible and thus the most interesting states, can be approximately described 
by wave functions like (12. 2p (the results which we are going to discuss in this section do not require 
strict factorization as in (|2.2p ; instead condensation of the initial wave function in the sense of (|3.ip 
would be sufficient). 

Because of the interaction among the particles, the factorization (|2.2p is not preserved by the 
time evolution; in other words, the evolved wave function ipN,t is not given by the product of one- 
particle wave functions, if t ^ 0. However, due to the mean-field character of the interaction each 
particle interacts very weakly (the strength of the interaction is of the order 1/N) with all other 
(N — 1) particles (at least in the initial state, every particle is described by the same one-particle 
orbital; every particles therefore "sees" all other particles). For this reason, we may expect that, in 
the limit of large N, the total interaction potential experienced by a typical particle in the system 
can be effectively replaced by an averaged, mean-field, potential, and therefore that factorization is 
approximately, and in an appropriate sense, preserved by the time evolution. In other words, we 
may expect that, in a sense to be made precise, 

N 

ip N)t (xi,...,x N )^l\_<Pt(xj) as N -> co (2.3) 

i=i 

for an evolved one-particle orbital (ft ■ Assuming (|2.3[) , it is simple to derive a self-consistent equation 
for the time-evolution of the one-particle orbital (ft- In fact, (j2.3f) states that, for every fixed time t, 
the N particles are independently distributed in space with density ^(a:)! 2 . If this is true, the total 
potential experienced, for example, by the first particle can be approximated by 



1 ^ f N — 1 

-J2v(x 1 -x j ) ~ -J2 / dyV( Xl -y)\My)\ 2 = ~ 1 ^( V * M*) * ( V 



*M 2 ) 



j>2 j>2 



as N —>■ cxd. It follows that, if (|2.3p holds true, the one-particle orbital (ft must satisfy the self- 
consistent equation 

id t ip t = -A(p t + (V*\tp t \ 2 )ipt (2.4) 



with initial data (ft=o = <*p given by (|2.2p . Equation (|2.4p is known as the Hartree equation; it is an 
example of a cubic nonlinear Schrodinger equation on M. d . Starting from the linear Schrodinger equa- 
tion (jl,2p on R rfAr , we obtain, for the evolution of factorized wave functions, a nonlinear Schrodinger 
equation on M. d ; the nonlinearity in the Hartree equation is a consequence of the many-body effects 
in the linear dynamics. 

The convergence of ipNt to the factorized wave function on the r.h.s. of (|2.3p as N — > cxd cannot 
hold in the L 2 -sense; we cannot expect, in other words, that \\ipN,t — <pf \\ —>■ as N — > oo (we 
use here the notation 93®^ (x) = Y[j=i f( x j))- Instead, (|2.3p has to be understood as convergence of 

(k) 

marginal densities. Recall that, for k = 1, . . . ,N, the A;-particle marginal ^ Nt associated with i^N,t 



is defined as the non-negative trace class operator on L 2 (M. dk ) with kernel given by 

7^ t ( x fc; x fc)= / dxN-kiN,t( x k,XN-k;x k > x N-k) 

I _ ^ (2.5) 

dxAf.fc -0Af,t(xfc, *N-k) V'AT.tCx'fe, *N-k) ■ 

It turns out that (|2.3p holds in the sense that, for every fixed k € N, the fc-particle marginal density 
associated with the left hand side converges, as N — > 00, to the fe-particle marginal density associated 
with the right hand side (which is independent of N, if N > k). In other words, assuming (|2.2p . one 
can show that, for a large class of interaction potentials, and for every fixed t € R and k 6 N, 



Tr 



7/fl " \<Pt)(<ft\ m ^0 as TV -► cxd , (2.6) 



where (ft is the solution to the Hartree equation (|2.4p with initial data 934=0 = </?• It is important 
here that the time t E M and the integer k > 1 are fixed; the convergence is not uniform in these two 
parameters. From (|1.5j) . we observe that (|2.6p implies (and is actually equivalent to the condition) 
that, for every fixed t£t and k £ N, and for every fixed compact operator J^ k > on L^( 



(lfar,t, («/ (fc) ® l (iV - fc) ) ^) - <vf*. ^ (fc Vf ) (2-7) 

as iV — > cxd. This means that, if we are interested in the expectation of observables which depend 
non-trivially on a fixed number k of particles, then we can approximate, as iV — ► 00, the true 
solution ij)N,t to the iV-body Schrodinger equation by the product of iV copies of the solution <pt to 
the Hartree equation (|2.4p . This approximation, however, is in general not valid if we are interested 
in the expectation of observables depending on a macroscopic (that is, proportional to N) number 
of particles. 

The first rigorous results establishing a relation between the many body Schrodinger evolution 
and the nonlinear Hartree dynamics were obtained by Hepp in [21] (for smooth interaction potentials) 
and then generalized by Ginibre and Velo to singular potentials in [20]. These works were inspired 
by techniques used in quantum field theory. We will discuss this method in Section UJ where we 
present a recent proof of (|2.6p . obtained in collaboration with I. Rodnianski in [28], which provides 
a quantitative control of the rate of convergence and makes use of the original idea of Hepp. 

The first proof of the convergence (|2.6p was obtained by Spohn in [30], for bounded potentials. 
The method introduced by Spohn was then extended to singular potentials. In [16], Erdos and Yau 



proved (|2.6[) for a Coulomb potential V(x) = ±l/|x|; partial results for the Coulomb potential were 
also obtained by Bardos, Golse and Mauser in [5] (note that recently a new proof of (|2.6p for the 
case of a Coulomb interaction has been proposed by Frohlich, Knowles, and Schwarz in [H5]). In [9|, 
a joint work with A. Elgart, we considered again the Coulomb potential, but this time assuming a 
relativistic dispersion for the bosons. Recently, in a series of papers [11 |, ll!2 |, ITH ], I14 (, I15j in collaboration 
with L. Erdos and H.-T. Yau (and also in [Sj, a collaboration with A. Elgart, L. Erdos and H.-T. 
Yau) the strategy of [30] was applied to systems with an A-dependent interaction potential, which 
converges, in the limit A — ► oo, to a delta- function. Note that in the one-dimensional setting, 
potentials converging to a delta-interaction have been considered by Adami, Golse and Teta in [2] 
(making use of previous results obtained by the same authors in collaboration with Bardos in pQ) . 
We will discuss these systems in Section [3j 

Recently, a different approach to the proof of (|2.6p has been proposed by Frohlich, Schwarz and 
Graffi in [17]. For smooth potentials, they can consider the mean-field limit uniformly in Planck's 
constant K (up to errors exponentially small in time); this allows them to combine the semiclassical 
limit and the mean field limit. It is also interesting to remark that the mean-field limit (|2.6j) can be 
interpreted as a Egorov-type theorem; this was observed by Frohlich, Knowles, and Pizzo in [18]. 

2.2 Derivation of the Hartree Equation for Bounded Potentials 

We consider, in this section, the dynamics generated by the mean field Hamiltonian (|2.ip under the 
assumption that the interaction potential is a bounded operator. We will assume, in other words, 
that V G L°°(R ) (recall that the operator norm of the multiplication operator V(xi — Xj) is given 
by the L°°-norm of the function V). To simplify a little bit the analysis we will also assume the 
external potential V ex t in the Hamiltonian (12. lj) to vanish; the techniques discussed here can however 
be easily extended to V ext / 0. 

Theorem 2.1 (Spohn, [30]). Suppose that 

N 1 

H N = J2 -Aj + - J2 V( Xi - Xj ) 

3=1 i<j 

with V G L°°(R d ). Let Vjv = <p 9N G L 2 (R dN ) for some tp G L 2 (R d ) with \\<p\\ = 1. Let tp Njt = 
e~ iHNt ij)N, and denote by j^ N ' t the k-particle marginal density associated with ipN,t- Then, for every 
fixed i£l, and for every fixed k > 1, we have 

Tr\^ t -\ip t )(^ k \ ^0 (2.8) 

as N — > oo. Here (ft denotes the solution to the Hartree equation 

id m = -Aip t + (y*\ip t \ 2 )p t (2.9) 

with initial data <ft=o = <*P- 

(k) 

Proof. The proof is based on the study of the time evolution of the marginal densities "f N t in the 
limit N — > oo. From (|1.2p . it is simple to show that the dynamics of the marginals is governed by a 



hierarchy of N coupled equation, commonly known as the BBGKY hierarchy: 

k , k 



^{S = E[-^7iS] + ^E[^( 






7=1 i<j 

k (2-10) 

(N - k) s-^ r ( fc+ i)- 

H ^7 2^ Tl ' fc + 1 V ( X J ~ X k+V^N,i ' 



N 

3=1 



We use here the convention that j Nt = if k > N. Moreover [A, B] = AB — BA denotes the 
commutator of the two operators A and B. The symbol Tr^+i denotes the partial trace over the 
(k + l)-th particle; the kernel of the /c-particle operator Tr^+i [V(xj — Xk + i),"/ Nt ] is given by 



Tr 



fc+i 



f2.ll) 



V(xj - x k+ i), 7SJ^ 1} J (xa.; x' fc ) 

da^fc+i (V"(wj - Xk+x) - V(x'j - Xfe+i)) 7 (fc+1) (x fc ,Xfc + i;x' fc ,x fc+ i) . 
Rewriting the BBGKY hierarchy (|2,10p in integral form, we find 
7$ = M«(t)7(*) + 1 jf d,W«(t _ fl)A W 7 W + (l - £) jf d. W«(< - s)^*^ (2.12) 

where U^ >{t) denotes the free evolution of A: particles, defined by 

z/ fc ) (t)7 « = e i*£? =1 A, 7 (fc) e -it£j =1 A, (2 13) 

and the maps A( fe ) and ZJW are defined by 






(fc) 



(2.14) 



and, respectively, by 



£(*) 7 (fc+D = -i^Tr fe+1 [y(x, - x fc+1 j, - 



(fc+i) 



(2.15) 



Note that B^ k > maps {k + l)-particle operators into fc-particle operators (while A^ k ' maps A:-particle 
operators into fe-particle operators). Since we are interested in the limit N — ► oo with fixed k > 1, it 
is clear that the second term on the r.h.s. of (|2.12p . as well as the contribution proportional to k/N 
to the third term on the r.h.s. of (|2.12p should be considered as small perturbations. Iterating the 
integral equation f)2. 12|) for n times, and stopping the iteration every time we hit a perturbation, we 
obtain the Duhamel type series 

n-l rt 

+ Wd Sl .. 







Sm— 1 



?Tl=l 

ft fSn-1 

+ / dsi . • • / 

Jo Jo 

1 N /•* 



771=1 







d Sm ^W(t - Sl )B^U^ +1 \ Sl - s 2 ) ■ • • J B( fc+m - 1 ^( A;+m )( Sm )7S fc+m) 

d Sn ^( fc )(t - sl )B«w( fe+i )( S i - S2 ) . . . s^-^i 



Sm, — 1 



^t + m-l^ ...p- dSmWW(t _ 8l)B (« 



771=1 



2V 



R (fc+m-lWfc+m) 



(2.16) 



To show (12. 8j) . we need to compare 7Jy t with 7 



(fc) 



00, t 



I^X^t 



where cpt is the solution to the 



Hartree equation (|2.9p . It is simple to check that the family {7^ t }k>i solves the infinite hierarchy 



(*) 



(written directly in integral form) 

which leads, after iteration, to the expansion 



(2.17) 



/o 

n-l rt 

l 

m=l JG /n 

ft j-- :: 

ds\ . , 



Sm— 1 



+ ^ / dsi... / 'd Sm ^)(t- Sl )^ fc )^ +1 )( Sl - S2 )...S( fc+m - 1 W( fc+m )( Sm ) 7 S fc+m) 

" ds n W (fc) (i - si)5Ww( fe+1 )(si - s 2 ) • • • 5( fc+n - 1 ) 7 ^+") 

(2.18) 



The difference between <y Nt and 7^ t can thus be bounded by 



Tr 



(ft) _ (fc) 
7jv,t 7oo,t 



ft fSn-l 

< / dsi . . . / 
Jo Jo 



ds n Tr 



«i*'(i - S i)b"=»w"= +i »(»i -«>... b'^"-" (t?,:; 1 - 7SS) 



4£jO- 



m=l 



ds m Tr 



+ ^ ^ / dsi... / d*„, Ti- 



m=l 



AT 



«W(t - si)BW . . . J B (fc+m ~ 1) 7^r ) 



Next we observe that, since U^ k '(i) is a unitary operator, 



Tr 



tf(*)(t) 7 (*) 



Tr 



7 



(*) 



Moreover, since V is a bounded potential, we have 



Tr 



< AT II VII Tr 



and 



where we used the fact that 



Tr 



S( fc ) 7 (*+ 1 ) < 2A;||F||Tr 



(ft) 



Tr 



Tr fc+ i7 



(fc+i; 



< Tr 



7 
(fc+i) 



(fc+i; 



(2.19) 
(2.20) 

(2.21) 

(2.22) 
(2.23) 



(Here the trace on the r.h.s. is a trace over (k + 1) particles.) Applying these bounds iteratively to 
the terms on the r.h.s. of (|2.19p . and using the a-priori information Tr 

analogously for 7^ t n> ), we obtain 



(fc+n) 

TN,t 



Tr 7 ^ n) = l(and 



Tr 



(ft) _ (ft) 
7iv,i 7oo,t 



m'(Jfe-l)! 



^ ' m=l 

Z,ofc+l ^ 

<2 fc (4||viitr + ^^£(4||viiir. 



(2.24) 



A^ 



m=l 



If < t < t , with t = 1/(8||V||), it follows that 



Tr 



(fc) _ (fc) 

7jV,t 'OQ,t 



< — h 

- 2« 



k 2 k+1 

N 



Since the l.h.s. is independent of the order n of the expansion, it follows that 

k2 k+l 



Tr 



(k) _ (k) 
^N,t 'oo,t 



< 



and thus that 



Tr 



(k) _ (k) 

7jv* 7oo,t 



for all < t < to an d for all k > 1. Next, set 



ti := sup < t > : lim Tr 



!N,s !oo,s 



N 
as N — ► oo 

for all fixed < s < t and fc > 1 



(2.25) 
(2.26) 



From (|2.26p . it follows that t\ > to- We show that t\ = oo by contradiction. Suppose that t\ < oo. 
Then, if t2 = ti — (to/2), we have, by definition, 



lim Tr 



(fc) _ (k) 
7AT,t 2 7oo,t 2 



for all k > 1. 



Starting from (|2.27f) . we are going to prove that 



lim Tr 

N— >oo 



(fc) _ (fc) 

7jv,t 7oo,t 



(2.27) 



(2.28) 



for all k > 1 and for all < t < ti + (to/2); this contradicts the definition of t\. To show (|2.28p . 
we expand ^y Nt and 7^ in Duhamel series similar to (12.16 
*i ~~ (*o/2)- Analogously to (|2.19p . we obtain, for t = t% + r, 

Tr|7^1-72i 
<Tr|^) ( r)( 7 g 2 - 7 S 2 ) 



we expand 7]^ t and 7^ in Duhamel series similar to (I2.16P and (I2.18p . but starting at time t2 



n-l 



+ E d ^ 



m=l 



ds r 



x Tr 



[T fSn-l 

+ / dsi... / ds n Tr 
Jo Jo 

+ -T7 E / dsi... / d.*,„'J.V 



m=l 



m=l 



N 



W«(t - Sl )i?( fe ^( fe+1 )( Sl -*,)... i?( fe + m - 1 ^( fc + m )( Sm ) ( 7 ^ 2 n) - 7 2£°) 
W«(t - Sl )^ (fc) • • • £ (fc+ - 1} (7^;° " 7&J?) 

W (*)( r _ ai ) 5 (*) . . .^( fe +— D( Sm _ 1 _ Sm )^( fc+m - i ) 7 ^r 1) 

W (*)( T _ Sl ) S W . . . ^(fc+ m -i) 7 Cfc+«x) 



'\ fc _|_ m _ ]_ /">" fSm-l 



(2.29) 



With respect to (12.19}) . we have one more term on the r.h.s. of the last equation, due to the fact that 
at time t = ti the densities do not coincide (while they do at time t = 0). Analogously to (|2.24p we 
find 



Tr 



(*) _ (*) < 2 k V — Tr 
m=0 



(fc+m) (fc+m) 

7jV,t 2 ~~ 7cxD,t 2 



+ 



fc2 fc+1 

N 
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if t\ — (to/2) < t < t\ + (to/2) (that is, if < r < to)- Choosing first n > sufficiently large (to make 
the second term on the r.h.s. smaller than e/3), and then N > sufficiently large (this guarantees 
that the third term, and, by (|2.27p . also the first term, are smaller than e/3), the quantity on the 
l.h.s. can be made smaller than any e > (for arbitrary k > 1 and t\ — (to/2) < t < t\ + (to/2)). 
This shows (|2.28|) and completes the proof of the theorem. □ 

2.3 Another Proof of Theorem 12.11 

From the proof of Theorem I2.1l presented above, we notice that the expansion of the BBGKY hierar- 
chy in (|2.16p is much more involved than the corresponding expansion (12.180 of the infinite hierarchy 
(|2.17p . It turns out that it is possible to avoid the expansion of the BBGKY hierarchy making use 
of a simple compactness argument; this will be especially important when dealing with singular po- 
tentials. In the following we explain the main steps of this alternative proof to Theorem 12.11 Then, 
in the next section, we will illustrate how to extend it to potentials with a Coulomb singularity. 

The idea, which was first presented in [5] HI [16] , consists in characterizing the limit of the densities 
7Jy t as the unique solution to the infinite hierarchy of equations (|2.17j) ; combined with the compact- 
ness, this information provides a proof of Theorem l2.1i More precisely, the proof is divided into three 
main steps. First of all, one shows the compactness of the sequence {^ N t ]k>i with respect to an 

appropriate weak topology. Then, one proves that an arbitrary limit point {7^ t }fc>i of the sequence 

(wtlfc=l * s a s °l u ti° n t° t ne infinite hierarchy (|2.17p (one proves, in other words, the convergence 
to the infinite hierarchy) . Finally, one shows the uniqueness of the solution to the infinite hierarchy 
(|2.17p . Since it is simple to verify that the factorized family {7^ t }fe>i, with 7^ = \(p t }(tp t \® k for all 

k > 1, is a solution to the infinite hierarchy, it follows immediately that 7^ — > \ipt){<£t\ as N — » 00 
(at first only in the weak topology with respect to which we have compactness; since the limit is an 
orthogonal rank one projection, it is however simple to check that weak convergence implies strong 
convergence, in the sense (12. 8p ). Next, we discuss these three main steps (compactness, convergence, 
and uniqueness) in some more details. 

Compactness: Let C\ = i2 1 (L 2 (IR dfc )) denote the space of trace class operators on L 2 (R dk ), 
equipped with the trace norm 

plli = Tr \A\ = It (A*A) 1/2 for all A € C\ . 

Moreover, let /C^ = )C(L 2 (M. dk )) be the space of compact operators on L 2 (M. dk ), equipped with the 
operator norm. Then C\ and Kk are Banach spaces and C\ = KL* k (see, for example, [27] [Theorem 
VI. 26]). By definition, the fe-particle marginal density ^ Nt is a non-negative operator in C\, with 

ll7Klli=Tr|7Sf,il=Tr 7 ^ = l 



sequence { i y^ t }N>k is compact with respect to the weak* topology of C\. 



for all N > k. For fixed t £ R and k > 1, it follows from the Banach- Alaouglu Theorem that the 

-1 

-A:' 

(k) 

Since we want to identify limit points of the sequence 7^ t as solutions to the system of integral 
equations (|2.17p . compactness for fixed t E M is not enough. To make sure that there are subsequences 

(k) 

of "y Nt which converge for all times in a certain interval, we use the fact that, since /C/% is separable, 
the weak* topology on the unit ball of C\ is metrizable. It is possible, in other words, to introduce 
a metric % on C\ such that a uniformly bounded sequence {A n }„ 6 pj G C\ converges to A e C\ as 
n — » 00 with respect to the weak* topology of C\ if and only if % (A n , A) — ► (see [2U] [Theorem 3.16], 
for the explicit construction of the metric %). For arbitrary T > let C([0,T],C k ) be the space 
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of functions of t E [0, T] with values in C\ which are continuous with respect to the metric rjk] on 
C([0, T],C\) we can define the metric 

%(7 (fe) (-),7 (fc) (-)):= sup i7fc(7 (fc) (<),7 (fc) (*))- (2-30) 

te[o,T] 

Finally, we denote by r pro( i the topology on the space © fc>1 ^([0, T], C\) given by the product of the 
topologies generated by the metrics % on C([0, T],/^,). 

The metric structure introduced on the space © fc>1 C([0, T], C\) allows us to invoke the Arzela- 

Ascoli Theorem to prove the compactness of the sequence Ys,t = {7jvt}fe=l' ^ e obtain the following 
proposition (for the detailed proof, see, for example, [131 Section 6]). 

Proposition 2.2. Fix an arbitrary T > 0. Then the sequence Tnj = (TzviK^l ^ ©fc>i ^([0'^1'^fc) 
is compact with respect to the product topology T pro d defined above. For any limit point T^t = 
{loot}k>l> loot * s symmetric w.r.t. permutations, non-negative and such that 

Tr^] t <l (2.31) 

for every k > 1. 

Remark. Convergence of T^t = {lNt}k=i ^° ^°°,t = {7oo tlfc>i with respect to the topology 
r prod is equivalent to the statement that, for every fixed k > 1, and for every fixed compact operator 
J (fc) elC k , 

Tr J (fc) (7$ " 7S) - (2.32) 



as iV — ► 00, uniformly in t for £ E [0,T]. Compactness of rjv,t with respect to the topology r pro d 
means therefore that for every sequence {Mj}j^ there exists a subsequence {Nj}j^ C {Mj}j G n 
and a limit point r^^ such that Tpf-t - * r^ in the sense f)2.32j) . 

Convergence: The second main step consists in characterizing the limit points of the (compact) 
sequence r_/v,t = {7JvtH>i as solutions to the infinite hierarchy of equations (|2.1Tj) . 



Proposition 2.3. Suppose that V E L°°(M. ) such that V(x) — ► as \x\ — ► 00. Assume moreover 
that Tooj = {7o | t }fc>i E 0fc>iC'([O,T],£^) is a Zirrai point of the sequence T N j = {lN,t}k=i with 
respect to the product topology T pro( i. Then 7^ = \<p)(ip\® and 



t\=^ k \t)tl+ fdsU^{t-s)B^ht: t 1] (2-33) 

Jo 



for all k > 1. HerelA^ >{t), and B^ ' are defined as in \2. 13\) and, respectively, in \2.15\) . 

Note that in Proposition 12.31 we assume the potential to vanish at infinity. This condition, which 
was not required in Section 12.21 is not essential but it simplifies the proof and it is also satisfied for 
the singular potentials (like the Coulomb potential) that we are going to study in the next sections. 

Proof. Passing to a subsequence we can assume that rV,* — ► ^oo,t as N — > 00, with respect to the 
product topology r pro d; this implies immediately that 7oo,o = |y)(y| • To prove (|2.33p . on the other 
hand, it is enough to show that for every fixed k > 1, and for every fixed j' ' from a dense subset 
of K k , 

Tr J^iSt = Trj( fc ^( fc )(t) 7 2 ) + f dsU^(t - s)Tr jW^W^+D . (2.34) 



'0 
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To demonstrate f|2.34j) . we start from the BBGKY hierarchy (|2.10p which leads to the relations 

/o 



TrjW 7 «=Trj( fc )^)( i)7 ^ + ^^/ ds Tr, J^ k \t - s) Vfo-x^g 



i=i 



,(*) 



+ 



JV-A; 



2V 



f t dsTrj( k ^ k \t-s)B^+ 1 K 
Jo 



(2.35) 



Since, by assumption, the l.h.s. and the first term on the r.h.s. of the last equation converge, as 
N — ► oo, to the l.h.s. and, respectively, to the first term on the r.h.s. of 112. 34j) (for every compact 
operator J' fc '), (|2.33p follows if we can prove that 



1 

N 



i=i 



Y, dsTrjWuM(t-s) V{xi- Xj )^ 



.,(*) 



(2.36) 



and that 



N-k 

N 



dsTrJ^U^(t-s)B^+ 1) ^ / dsTrJ^U^(t - s)B^g+^ 



(2.37) 



as N — ► oo. Eq. (|2.36p follows because 

Trj(%W(t-s) [^-x^T^J <2||j( fc )||||F||Tr 



(fc) 



<2||J (fc) ||||F| 



is finite, uniformly in N. To prove Eq. (|2.37|) one can use a similar argument, combined with the 
observation that 

Trj(%W(t- S )^)( 7 ^- 7 ^)) 

= fcTr [(U( k \s - t)J«) V( Xl - x k+1 ) - V( Xl - x k+l ) (lA^{s - t) J«)] (ffi> - -ygj 1 ) 



as N — > oo. This does not follow directly from the assumption that Tjf,t ~> Too,t with respect to 
the topology r pro d because the operators (U^ k '(s — t)J^ ')V(x\ — Xfc+i) and V{x\ — Xk+i)(U( k >(s — 
t)j( k >) are not compact on L 2 (M. d ( k+1 >). Instead we have to apply an approximation argument, 
cutting off high momenta in the x^+i-variable, and using the fact that, by energy conservation, 
TrV^ +l7 ] vt Vfc + i is bounded, uniformly in N and in t (and that, therefore, TrV£ + i 7oo t *fe+i I s 
bounded as well). Note that, because of the assumption that V(x) —> as \x\ — > oo, we only need a 
cutoff in momentum, and no cutoff in position space is necessary. The details of this approximation 
argument can be found, for example, in Eq. (7.35) and Eq. (7.36) in the proof of Theorem 7.1 in 
[13] (after replacing 5p through the bounded potential V). □ 



Uniqueness: to conclude the proof of Theorem 12.14 we still have to prove the uniqueness of the 
solution to the infinite hierarchy (I2.33p . 

Proposition 2.4. Fix r^o = {7oooi fc>1 ^ ©fc>l^jfc- Then there exists at most one solution 



oo.t 



„( fc ) 



„( fc ) 



{Too tlfc>l ^ ©fe>i C{^,T\,C\) to the infinite hierarchy h2.33\) such that 7oo t=0 = Too o an( ^ 



(k) 



M 



Tr \loo t\ < 1 f° r allk>l and all t £ [0,7]. 
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Proof. Suppose that {7^ \ t }k>i and {7^ 2 t}fc>i are two solutions of (J2.33P with the same initial 
data {7ooo}fc>i' sucn that Tr^^J < 1, for all k > 1, t E [0, T], and for i = 1,2. Then we can 
expand 7^ \ t and 7^ 2 1 in the Duhamel series (|2.19p . It follows that 



Tr 



(k) _ (k) 

7oo,l,t 7oo,2,t 



< / dsi 

'0 



ds n Tr 



u^xt - s l)B ^ . . . flOH-D ( 7 ^ n - 7 2ri 



Applying recursively the bounds (|2.20[) and (|2.22j) . we obtain 

Tr 
and thus, for < t < l/8||y 



(k) _ (k) 

7oo,l,t 7oo,2,i 



-l^!^r (2||y||t)n - 2fe(4||y||tr 



Tr 



(k) _ (k) 

7oo,l,t 7oo,2,t 



<2 



k—n 



Since the l.h.s. is independent of n > 1, it has to vanish. This proves uniqueness for short time. 
Iterating the same argument, we obtain uniqueness for all times. □ 

2.4 Derivation of the Hartree Equation for a Coulomb Potential 

The arguments presented in Section T2.2I and in Section T2.3I required the interaction potential V to be 
bounded. Unfortunately, several systems of physical interest are described by unbounded potential. 
For example, in a non-relativistic approximation, a system of gravitating bosons (a boson star) can 
be described by the Hamiltonian 



H 



N 



N 



A 



N 

E 



1 






(2.38) 



with a singular Coulomb interaction among the particles. The factor of \/N in front of the potential 
energy can be justified, when describing gravitating particles, by the smallness of the gravitational 
constant. As in the case of bounded potential, we are interested in the dynamics generated by the 
Hamiltonian (j2.38[) on factorized initial A^-particle wave functions. We specialize here in the phys- 
ically most interesting case of particles moving in three dimensions; however, the theorem remains 
valid in all dimensions d > 2. 



Theorem 2.5 (Erdos-Yau, [S]). Let ipN 



Y 



®N 



for some (p E H 1 



and let ipjy-i 



-iHfft 



■IpN 



where the Hamiltonian Hn is defined as in 12.381) . Then, for arbitrary k > 1 and t EM, we have 



Tr 



(fc) 
TNt 



\<Pt)(<Pt\ 9k 



(2.39) 



as N — y 00. Here tpt is the solution to the nonlinear Hartree equation 

id t ipt = -Atp t - A ( — * \if t \ 2 ) (ft 



with initial data <ft=o = f- 

Remark. Although, physically, the value of the constant A is positive (corresponding to the 
Coulomb attraction among gravitating particles), the theorem remains valid also for negative values 
of A (corresponding to repulsive Coulomb interaction). 
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The general strategy used in [16] to prove Theorem 12.51 is the same as the one outlined in 
Section [2.31 First one proves the compactness of the sequence of marginal {TwilfeLi with respect to 
an appropriate weak topology (the product topology r pro d introduced after (|2.3U|) ). then one shows 

(k) (k) at 

that an arbitrary limit point {7^ t }k>i of the sequence {7]yt}fc=i ^ s a solution to the infinite hierarchy 
of equations 



7 (*) = W (fc)( 4 ) 7 (k) + r ds^ k \t ~ S )BW 7 f +i ) 

Jo 



(2.40) 



where US ' is the free evolution defined in (|2.13p . and the collision map B^ ' is now given by 



B (k) Ak+l) 






„(*+i) 



Xk+l\ 



(2.41) 



Finally, one proves the uniqueness of the solution to (|2.40p . Although the proof of the compactness 
and of the convergence also require several changes with respect to what we discussed in Section [273], 
the main difficulty one has to face when the bounded potential is replaced by the Coulomb interaction 
is the proof of the uniqueness of the solution to the infinite hierarchy. The key idea introduced by 
Erdos and Yau in |16| was to restrict the class of densities for which uniqueness must be proven. In 
Proposition 12.41 uniqueness is proved in the class of densities with Tr|7 t | < 1 for all k > 1, and 
all t S [0,T] (but the same argument works under the weaker assumption TY^ | < C , for some 
constant C < 00). Following [16], in the case of a Coulomb potential we are only going to show the 
uniqueness of (|2.40p in the class of densities T t = J7 t }k>i satisfying the a-priori bound 



Tr 



(1 - A^ 2 ... (1 - A fc ) 1/2 7 f } (1 - A fc )V2 . . . (i _ Al) i/2 



<C k 



(2.42) 



for all k > 1 and for all t G [0, T\. Note that, for non-negative densities 7^ > (in the sense of 
operators, that is, in the sense that {if)( k \^l 'ip^) > for all ip^ G L 2 (R 3h )) we have 



Tr 



(1 - A1) 1 / 2 ... (1 - A,) 1 /^) (1 - A fc )V2 . . . (i _ Al ) 1/2 = Tr (1 - A0 • . . (1 - A fc ) 7t ( 



» 



There is, of course, a price to pay in order to restrict the proof of the uniqueness to this class of 
densities. In fact, to apply this uniqueness result to the proof of Theorem 12.51 one h as to show that 



an arbitrary limit point T 



oo,t 



„(*) 



W-lJV 



{7oo t\k>i of the sequence of densities Tnj = {7_/vt}fc=i associated 
with 1^1 N,t satisfies the a-priori bound (|2.42p . Due to the Coulomb singularity, this is actually not so 
simple and requires an additional approximation argument. 

Approximation of the Coulomb singularity: For a fixed e > we define the regularized 
Hamiltonian 



H 



N 



N 

E 



A N 

-Y 



1 









Moreover, for a fixed sufficiently small 5 > 0, we introduce the regularized initial data 

X{5H N /N)^ N 



f 



N 



\ X (6H N /N)^ N \ 



(recall that tp^ = ip 



®N\ 



(2.43) 



(2.44) 



where x G Co°(^) i s a monotone decreasing function such that x( s ) = 1 f° r all s < 1 and x( s ) = 
for all s > 2. We consider then the regularized evolution of the regularized initial wave function 



f/jN,t 



-iH^ft 



i> 



N ■ 
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The advantage of working with the regularized wave function tpN,t instead of ipN,t is that it satisfies 
the following strong a-priori bounds. 

Proposition 2.6. Let ipiyj = e~ tHNt i^i^, for some fixed e, 5 > 0. Then there exists a constant C > 
(depending on e, 5) and, for all k > 1, there exists Nq = No(k) > k such that 



for all N> N Q . 



(^,t,(l-A 1 )...(l-A fc )Viv, i )<C A 



#) 



(2.45) 



Remark. Expressed in terms of the /c-particle marginal ^y Nt associated with tpNj, the bound 
(g3SD reads 

Tr (1 - Ax) ... (1 - A fc ) 7^1 < C k . (2.46) 

We will show Proposition 12.61 below, making use of Proposition 12.7} we will see there that the 
regularization of the Coulomb singularity and of the initial wave function both play an important 
role. For the solution ipN,t of the original Schrodinger equation with the original factorized initial 
data tpN = ip® N it is not known whether bounds like (|2.45p hold true. 

In order for the regularized wave function ipN,t to be useful, one needs to prove that it approxi- 
mates, in an appropriate sense, the wave function ipN,t- To compare the two TV-particle wave function, 
we introduce a third wave function i^N t = e~ lHNt tpN, and we use the triangle inequality 



•>N. 



t ~ ll>N,t\\ ^ WN,t ~ $N,t\\ + WN,t ~ lpN,t\\ ■ 



(2.47) 



The second term is actually independent of time because of the unitarity of the evolution. Using the 
definition of the regularized initial data ipN, one can prove that 

life - Mt\\ = Un - $n\\ < C6 1 / 2 
uniformly in N. To control the first term on the r.h.s. of (|2.47p . we observe that 







d 
dt 


1pN,t ~ $N,t 


2 


= 2bn ((h n - H^ $ Nit 


, 1pN,t ~ 1pN,t 


} 


and thus that 




d 
ch 


lpN,t ~ 4>N,t 


2 


< 2 


[H N - H N J \j) N ,t 




1pN,t ~ i>N,t 




We have 


110? 


In- 


H N )^N,t\\ = 




AT2Z 
i 


^ 1 
> 1 |7j ; 

\Xi Xj ( \Xi C 


"j 


+ eN-t)'*"' 


'.t 



(2.48) 



Using the permutation symmetry of ipN,t, it follows that 



(H N -H N )^ Nyt \\ 
< £ 2 (V>jv,t, 



1 



1 



x\ — x 2 \ i\x\ - x 2 \ + eN x ) \x 3 - x 4 | (|x 3 - x 4 | + eN 

1 



ZitV'JV,, 



+ e 2 N- l U NyU -j 

\ \X\ - X 2 \ (|Xl - 



x 2 \ + eN- 1 ) \x 2 - x 3 \ (\x2 - x 3 \ + eN x ) 



riT W> ; 



+ e 2 N- 2 (i; N t, 



1 



xi - x 2 \ 2 (\x\ - x 2 \ + eN x ) 
16 



tpN;. 



and thus 



(H N -H N )il> Nj tf<e 2 U N 



\xi — X 2 \ 2 \x3 — xA 2 



t 



N.t 



+ e 2 N- 1 ($M t 



' \x± - x 2 \ 2 \x 2 - x 3 \ 2 



1>. 



N,t 



+ e l l 2 N^I 2 U N „- 1 —^Jn,). 

\ \X\ — X 2 \ ' I 



Applying Hardy inequalities in the form 

1 



X 7 X o 



^Cil-Aif^il-Aj^/ 2 



(2.49) 



for every < a < 3, if /3 + 7 > a (see [9] [Lemma 9.1] for a proof) we find that 

\\{H N - H N )$ N ,t\\ 2 < e 2 (l + N- 1 )^, (1 - Ai)(l - A 3 )fet) 



+ £ 1 /2 iV -l/2 ( ^ t) (l_ Al )( 1 



jM 



2)WN,t 



and thus, from ()2.46p . 



uniformly in iV". From (|2.48j) . applying Gronwall's Lemma, it follows that 

life- fell <cfe 1/4 t- 



From (|2.47p . we obtain that 



and thus 



Tr 



N,t-^NA\<c(e l lH + 5 l l 2 ) 
<c ( £ i/4 t + S i/2 



(fe) ~{k) 

Twt ~ r N t 



(2.50) 



for every k G N, uniformly in iV > fc. Because of ()2.50p . it suffices to prove (J2.39P with ") Nt (the 
A;-particle marginal associated with ipN,t) replaced by 7^ (the /c-particle marginal associated with 
the regularized wave function ipN,t) for fixed e, 5 > 0; at the end (|2.39p follows by letting e, 5 — > 0. 
Note that in [16] a slightly different approximation of the initial data was used; the details of the 
approximation presented above can be found (for a different model) in [12] [Section 5]. 

Energy estimates: To prove the a-priori bounds of Proposition l2.6|, one can use so called energy 
estimates; these are estimates that compare the expectation of high powers of the Hamiltonian with 
corresponding powers of the kinetic energy. 

Proposition 2.7. Suppose that Hn is defined as in {2.1^3^ with A > (the case A < is simpler). 
Then there exist constants C\ > 1 and C 2 > and, for every k > 1, there exists an Nq = No(k) G N 
such that 

& N , (H n + dN) k ,I> N ) > C%N k (Vv, (-Ax + Ci) . . . (-A fe + d)il> N ) (2.51) 

for every rpN G L 2 (M. ) (symmetric with respect to permutations). 
Proof. Using the operator inequality 

1 7T, , 

1 1 < T v i 

X. — T -\ A. 

1 " J 1 \ ■*■ 
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we can find a constant C\ > 1 (depending on the coupling constant A) such that 

]— ^ < \ (~^ + Ox) = Ul (2.52) 



where we defined Sj = (— Aj + Ci) 1 ' 2 . Note also that, for every < a < 3 there exists a constant 
C Q < oo such that 

I l -^<C a S«. (2.53) 

We are going to prove (|2.5ip for C\ fixed as in (|2.52|) . and for an arbitrary < C 2 < 1/2. The 
proof is by a two-step induction over k > 0. For k = 0, the claim is trivial. For k = 1, it follows from 
(|2.52p because, as an operator inequality on the permutation symmetric space L 2 (\R. 3N ), we have 

(H N + dN) > NSJ - — -. r > C 2 NSl . (2.54) 

2 \X\ — X 2 | 

Next we assume that (|2.5ip holds for all k < n and we prove it for k = n + 2, for an arbitrary n G N. 
To this end, we observe that, because of the induction assumption, 

(H N + CiiV) n+2 = (H N + C X N){H N + C 1 N) n (H N + dN) 
> C%N n (H N + CiA^)5 2 . . . Sl{H N + CiiV), 

for all iV > iV (n). Writing 

A ^ 1 

(^+djv)= E ^+ fc ^ with ^ = E^-^E b ._ x .| + £jV -i 

j>n+l j=l i<j l l Jl 



it follows that 



(H N + dN) n+2 > C%N n (N - n)(N -n- l)Sf . . . S 2 n 



+2 

4c2 c2 



+ C%N n {N - n)SfS 2 . . . S 2 n+1 (2.55) 

+ C^N n {N - n) (Sf . . . S 2 n+1 h N + h.c.) . 

The first two terms are positive. As for the third term, by the definition of /tjv, we find that 

(<fl q2 , ,, w (N-n)(N-n-l) ( 2 2 A 

(S 1 ...S n+1 h N + h.c.)> [ S i- S ^ \ Xn+2 - Xn+3 \ + eN-i +h - C - 

(N - n)n / 2 2 A , hc 

t>i • • ■ J n 4-i 1 ; : ~ ~ t ti-C. 1 . 

(2.56) 



2iV V n+1 |xi-x 2 |+eW- 1 



The first term on the r.h.s. of (|2.56|) can be bounded by 



Si ... S n+ x -. j— ■ — -— j- + h.c. < 2S\ . . . S n+ \- r 5,1+1 ■ ■ ■ S\ 

\x n+2 -x n+3 \ +eN x J |x n + 2 -x n + 3 | 



(2.57) 



-i D l • • • D n 



+2 
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As for the second term on the r.h.s. of (|2,56p . we remark that 
S " ■ ■ ■ S " +1 \ Xl - Xn+2 \ +£ N-i + h ' C 



= S n +i ■■■S 2 ( (-Ai + Ci), — — -— r + h.c. ) S 2 ■ ■ ■ S n+1 

V \x\ -x n+2 \ +eN l ) 

< 2Cl S n+ l . . . S 2 - — ——rS2 ■ ■ ■ S n+ l 

\xi - x n+2 \ +eN 1 

+ 25 n +l ... 52 Vi, j— ; — rVi^ . . . S n+ i 

\xi - x n+2 \ +eN L 



(xi - x n+2 ) 

\xi - x n+2 \ (\xi - x n+2 \ + eiV- 1 ) 



+ A 5 n+ i ...S 2 I V*- — ■ o + h.c. S 2 . . . S n+ i 



Applying a Schwarz inequality in the last term, we conclude that there exists a constant D > 
(depending on A) such that 

( 5 l • • • 5 '+i i T I \+ fN -i + h - c <DSf... S 2 n+2 . (2.58) 

V \X\ — X n +2| + £l ^ ' 

Similarly, using the operator inequalities (|2.53p . the last term on the r.h.s. of (|2.56j) can be bounded 
by 

(S\ . . . S 2 n+1 _ | + h.c.) < De- l N Sf... S 2 n+1 + DSfS 2 2 . . . S 2 n+1 (2.59) 

for all < e < 1 and for a constant D depending only on A (it is at this point that the condition 
e > is needed). Inserting (|2.57|) . (|2.58|) . and (|2.59|) in the r.h.s. of ([2.56J) . and the resulting bound 
in the r.h.s. of (|2.55p . we obtain that there exists Nq > (depending on n) such that 

(H N + dN)^ 2 > C% +2 S 2 . . . S 2 n+2 

for all A^ > Ao- Note that the value of Nq also depends on the parameter e > 0. □ 

Using the result of Proposition 12.71 it is simple to complete the proof of the a-priori bounds for 
Vw,* = e~ lHNt ipN (recall the definition of the regularized initial data tp^ hi (|2.44p ). 

Proof of Proposition HUH . From (I2.51J) . and since C\ > 1, we have 

$ N)t , (1 - Ai) . . . (1 - A fc )fet> < $N,t, (Ci - Ax) . . . (d - A k )$N,t) 

<-^k$N,t>(HN + CiN) k i; N , t ) 

' $ N ,(H N + CiN) k $ N ) 



C$N k 



where in the last line we used the fact that the expectation of any power of Hn is preserved by the 
time-evolution. From the definition ([2.44J) of t/'jVj we immediately obtain f|2.45[) . D 

Since the a-priori bounds for ^ N t obtained in Proposition 12.61 hold uniformly in N, they can also 
be used to derive a-priori bounds on the limit points {7^ t }k>i of the sequence { , y Nt }^ =1 - 
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„( fc ) 



Corollary 2.8. Suppose thatF^^ = {joot}k>l G 0fc>iC([O)^]i^) ^s a limit point of the sequence 



&)\N 



S k ) 



^N,t = {7jvt}fe=i w ^ respect to the product topology T pro d defined after \2.30\) . Then %J t > and 
there exists a constant C such that 



Tr(l-A 1 )...(l-A k ) 1 ^<C k 



(2.60) 



for all k > 1. 



Uniqueness: The bounds of Corollary 12.81 are crucial; from (|2.60j) it follows that it is enough to 
show the uniqueness of the infinite hierarchy (|2.33[) in the class of densities satisfying (|2.60p , a much 
simpler task than proving uniqueness for all densities with Tr |7 4 | < C k . 

Theorem 2.9. Fix {7^}fc>i G ® k >i£\- Then there exists at most one solution {7 t }fc>i G 
fe>1 C([0, T],£^) to the infinite hierarchy \2.40\) , such that 



Tr 



1 - A^ 2 ... (1 - A,)^ 7 W(l - A fc ) 1/2 . . . (i _ Al) i 



/2 



<c k 



(2.61) 



/or all k > 1, and aZ/ i G [0, T]. 
Proof. We define the norm 



7 (fc) llw fc =Tr 



(1 - A X ) 1/2 ... (1 - A^/y^l _ Afc ) 1/2 . . . (i _ Al ) X 



/2 



and we observe that there exists a constant C > with (recall the definition (12.411) for the collision 
map B^) 

\\B^H k+1) \\n k <CkU k ^\\ Hk+1 . (2.62) 

To prove (|2.62p . we write 

\\B^ k H k+1) \\H k < X> S 1 ...S k (Tr k+1] l - - 1 (^)s k ...S 1 



k 



Si . . . S k [ Tr fc+ i 7 



A k + l ) 



1 



\Xj Xfc_|_i| 



S k ... Si 



All terms can be handled similarly. We show how to bound the summand with j = 1 on the first 
line. 



Tr 



Si . . . S k \ Tr fc+1 



1 



\xi -x k+ i\ 



n 



(fc+i) 



s k ■ ■ ■ Si 



= Tr 
< Tr 



5*1 ... 5^ I Trfe-nS. , i 1 r 5,, 1 5fc + i7^ ^Sfc+i ) S k . . . Si 

V K+L \xi -x k+ i\ K+L J 



•Si^fc+li 



1 



- Ilsls *+ 1 ]^ : ^[ t ' / - 1 

— ^ II I lln-fc+i 



^k+i^i Si...S k S k +iT S k+ iS k ...Si 



9- 1 'T 1 !! IK( fc+1 )| 



n k 



+1 
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where in the second line we used the cyclicity of the partial trace, in the third line we used (|2.23[) 
and, in the last line, we used the bound 

fc+1 |a;i-a; fe+ i| fe+i 
To prove (J2.63P we write, assuming for example that k = 1, 

JlOo -. T<Sa 61 = <Si <Sn (1 — AlJl "An <Si 

Fi - x 2 \ \xi - x 2 \ 

5—1 Q — 1 Q — 1 Q— 1 

|a;i-x 2 | 

+ sr^v?: — - — rViS^sr 1 + s^s^vi , (xi ~ ^ s^sz 1 , 

Fl — ^2 1 \ x l — %2\ 

and we use the norm-estimates llViS 1 ^ || < 00 and ||S^~ |xi — a^l - "^ || < 00 for all < a < 2 (by 

Suppose now that {7- 1 }fc>i, for i = 1, 2 are two solutions to the infinite hierarchy (|2,40p . Using 
(J2.18D . we can expand both 7 J t and j 2 t m a Duhamel series. From (I2.62p . and from the fact that 

ll^ (fc) 7 (fc) llH fc = ll7 (fc) l|H fc , we obtain that 

(k) (k) (fc + n)! /•* r n_1 ^, iu(*+«) _ J*+») 



< ^^^ / <i*i . . . / d,„ h:;:> - 7 - „ - IIWfc+n 

JO JO 

< C fc (Ct) n 



for any n. Here we used the a-priori bounds (|2.6ip . For t < 1/(2C), the l.h.s. must vanish. This 
shows uniqueness for short time, and thus, by iteration, for all times. □ 

3 Dynamics of Bose-Einstein Condensates: the Gross-Pitaevskii 
Equation 

Dilute Bose gases at very low temperature are characterized by the macroscopic occupancy of a 
single one-particle state; a non-vanishing fraction of the total number of particles N is described 
by the same one-particle orbital. Although this phenomenon, known as Bose-Einstein condensation, 
has been predicted in the early days of quantum mechanics, the first experimental evidence for its 
existence was only obtained in 1995, in experiments performed by groups led by Cornell and Wieman 
at the University of Colorado at Boulder and by Ketterle at MIT (see O [6]). In these important 
experiments, atomic gases were initially trapped by magnetic fields and cooled down at very low 
temperatures. Then the magnetic traps were switched off and the consequent time evolution of the 
gas was observed; for sufficiently small temperatures, it was observed that the gas coherently moves 
as a single particle, a clear sign for the existence of condensation. 

To describe these experiments from a theoretical point of view, we have, first of all, to give a 
precise definition of Bose-Einstein condensation. It is simple to understand the meaning of conden- 
sation if one considers factorized wave functions, given by the (symmetrization of the) product of 
one-particle orbitals. In this case, to decide whether we have condensation, we only have to count the 
number of particles occupying every orbital; if there is a single orbital with macroscopic occupancy 
the wave function exhibits Bose-Einstein condensation, otherwise it does not. In particular, wave 
functions of the form ^)jv(x) = ]^[ =1 tp(xj), for some tp S L 2 (M 3 ) (we consider in this section three 
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dimensional systems only), exhibit Bose-Einstein condensation; since in these examples all particles 
occupy the same one-particle orbital, we say that ipw exhibits complete Bose-Einstein condensation 
in the state (p. 

Although factorized wave functions were used as initial data in Theorem 12.11 and Theorem 12.51 
they are, from a physical point of view, not very satisfactory, because they do not allow for any 
correlation among the particles. Since we would like to consider systems of interacting particles, 
the complete absence of correlations is not a realistic assumption. For this reason, we want to give 
a definition of Bose-Einstein condensation, in particular of complete Bose-Einstein condensation, 
that applies also to wave functions which are not factorized. To this end, we will make use of the 
one-particle density 7^ associated with an iV-particle wave function ip^. By definition (see (|1.4p ). 
the one-particle density is a non-negative trace class operator on L 2 (R 3 ) with trace equal to one. 



It is simple to verify that the eigenvalues of 7^ (which are all non-negative and sum up to one) 
can be interpreted as probabilities for finding particles in the state described by the corresponding 
eigenvector (a one-particle orbital). This observation justifies the following definition of Bose-Einstein 
condensation. We will say that a sequence {i/)n}n&n with ipN 6 L^(R 3iV ) exhibits complete Bose- 
Einstein condensation in the one-particle state with orbital (p € L 2 (M 3 ) if 

Tr | 7 «- 1^1 1 -0 (3.1) 

as N — > 00. In particular, complete Bose-Einstein condensation implies that the largest eigenvalue 
of 7^ converges to one, as N — > 00. More generally, we say that a sequence {iPn}ngn exhibits (not 
necessarily complete) Bose-Einstein condensation if the largest eigenvalue of 7^ remains strictly 
positive in the limit N — > 00. Note that condensation is not a property of a single iV-particle wave 
function tp^, but it is a property characterizing a sequence {^atj-jvcn in the limit N — > 00. 

It is in general very difficult to verify that Bose-Einstein condensation occurs in physically inter- 
esting wave functions of interacting systems. There exists, however, a class of interacting systems 
for which complete condensation of the ground state has been recently established. 

In [21], Lieb, Yngvason, and Seiringer considered a trapped Bose gas consisting of N three- 
dimensional particles described by the Hamiltonian 

N N 

H T P = Y, (" A i + V "*( x j)) + £ V ^ x i - Xj), (3.2) 

j=l i<j 

where V ext is an external confining potential with limui^oo V ex t(x) = 00, and V/v( x ) = N 2 V(Nx), 
where V is pointwise positive, spherically symmetric, and rapidly decaying (for simplicity, V can 
be thought of as being compactly supported). Note that the potential V/v scales with iV so that 
its scattering length is of the order 1/N (Gross-Pitaevskii scaling). The scattering length of V is a 
physical quantity measuring the effective range of the potential; two particles interacting through V 
see each others, when they are far apart, as hard spheres with radius given by the scattering length 
of V. More precisely, if / denotes the spherical symmetric solution to the zero-energy scattering 
equation 

-A + -V(x) I / = with boundary condition f(x) — ► 1 as |x| —>■ 00, (3.3) 



the scattering length of V is defined by 



ao = lim I re I — |x|/(x) . 

\x\— >oo 
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This limit can be proven to exist if V decays sufficiently fast at infinity. Another equivalent charac- 
terization of the scattering length is given by 

87ra = f dxV(x)f(x). (3.4) 

It is simple to verify that, if / solves (|3.3p . the rescaled function Jn(x) = f(Nx) solves the zero 
energy scattering equation with rescaled potential Vn, that is 

-A + -Vn ) In = with /n(x) —*■ 1 as \x\ — > oo . (3-5) 

This implies immediately that the scattering length of Vn is given by a = ao/N, where ao is the 
scattering length of the unsealed potential V. Note that, for \x\ 3> a, Jn(x) — 1 — a/|x|. For \x\ < a, 
Jn remains bounded; for practical purposes, we can think of this function as /jv(x) ~ 1 — a/(|x| + a). 

Letting N — > oo, Lieb, Yngvason, and Seiringer showed that the ground state energy E(N) of 
(|3.2[) divided by the number of particle N converges to 

E ( N ) „ / s 

hm — — — = mm £gp(<p) 

N^oo N v eL 2 (R3): 11^11=1 

where £qp is the Gross-Pitaevskii energy functional 

£ G p(<p)= [dx (|V^(x)| 2 + Kxt(x)|^(x)| 2 + 4vrao|^(x)| 4 ) . (3.6) 



Later, in [22J, Lieb and Seiringer also proved that the ground state of the Hamiltonian (j3.2|) 
exhibits complete Bose-Einstein condensation into the minimizer of the Gross-Pitaevskii energy 
functional <?gp- More precisely they showed that, if ipN is the ground state wave function of the 
Hamiltonian (|3.2p and if 7^ denotes the corresponding one-particle marginal, then 

7at -► I0Gp)(0Gp| asiV^oo, (3.7) 

where 4>qp £ L 2 (IR 3 ) is the minimizer of the Gross-Pitaevskii energy functional (j3.6h . 

To describe the experiments mentioned above, it is important to understand the time-evolution of 
the Bose-Einstein condensate after removing the external traps. We define therefore the translation 

invariant Hamiltonian 

N N 

j=l i<j 

and we consider solutions to the iV-particle Schrodinger equation 

idti> N ,t = H N i> N ,t => ^N,t = e- iHNt ^ N (3.9) 

with initial data ipN exhibiting complete Bose-Einstein condensation. In a series of joint articles with 
L. Erdos and H.-T. Yau, see \\.'2\ [T3| [Tl"} I15| . we prove that, for every fixed time t € R, the evolved 
A^-particle wave function ipN,t still exhibits complete Bose-Einstein condensation. Moreover we show 
that the time evolution of the condensate wave function evolves according to the one-particle time- 
dependent Gross-Pitaevskii equation associated with the energy functional £qp- Our main result is 
the following theorem. 
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Theorem 3.1. Suppose that V > is spherically symmetric and V(x) < C{x)~ a , for some a > 5, 
and for all x G M. 3 . Assume that the family {ipn}n&n with tp^ G L^(]R 3Af ) and \\iPn\\ = 1 f or a ^ N, 
has finite energy per particle, that is 

(4> N ,H N i> N )<CN (3.10) 

for all N € N, and that it exhibits complete Bose-Einstein condensation in the sense that 



Tr 



tS'-WM -0 (3.11) 



as N — ► oo for some ip G L (R ). Then, for every k > 1 and t£f, we /iaue 



TV 



,(*) 



7S-I^><^r -0 (3.12) 

as iV — > oo. ifere ^ zs t/ie solution of the nonlinear Gross-Pitaevskii equation 

id t <pt = -A<pt + 87ra |v?t| 2 ^t (3.13) 

TOi/i initial data <ft=o = <*P- 

Making use of an approximation of the initial iV-particle wave function (similarly to (|2.44p ). it 
is possible to replace the assumption (|3.10p with the much more stringent condition 

{ip N ,H k N ip N )<C k N k (3.14) 

for all k G N. In the following we will illustrate the main ideas involved in the proof of Theorem 13. 11 
assuming the initial wave function ipjy to satisfy (|3.14p . 

3.1 Comparison with Mean-Field Systems 



The formal relation of the Hamiltonian (|3.8p with the mean-field Hamiltonian (|2.ip considered in 
Section [2] is evident; (13. 8p can in fact be rewritten as 

N 1 

jr' = l i<j 

with vn(x) = N s V(Nx). Since we are considering three dimensional systems, vjy converges to a 
delta-function in the limit of large N; at least formally, as iV — > oo, we have vn(x) —> bo5(x), 
where bo = J V(x)dx. In other words, the Hamiltonian (|3.8p in the Gross-Pitaevskii scaling can be 
formally interpreted as a mean field Hamiltonian with an ./V-dependent potential which converges, 
as N — > oo, to a 5- function. Despite the formal similarity, it should be stressed that the physics 
described by the Gross-Pitaevskii Hamiltonian is completely different from the physics described 
by the mean field Hamiltonian (|2.ip . In a mean-field system, each particle typically interacts with 
all other particles through a very weak potential. The Gross-Pitaevskii Hamiltonian (J3.8I) . on the 
other hand, describes a very dilute gas, where interactions are very rare and at the same time very 
strong. Although the physics described by (|2.ip and (13. 8p are completely different, due to the formal 
similarity of the two models, we may try to apply the strategy discussed in Section 12.31 to prove 
Theorem 13.11 In other words, we may try to prove Theorem 13.11 by showing the compactness of 
the sequence r^t = {ijvtlfcLi w ith respect to an appropriate weak topology (it is going to be the 
same topology introduced in Section [2. 3p . the convergence to an infinite hierarchy similar to (|2.33p . 
and the uniqueness of the solution to the infinite hierarchy. It turns out that it is indeed possible 

24 



to extend the general strategy introduced in Section [2.31 to prove Theorem l3.lt however, as we will 
see, many important modifications of the arguments used for bounded or for Coulomb potential are 
required. We discuss next the main changes. 

First of all, the simple observation that, formally, vn(x) — ► bo5(x) as N — ► oo may lead to the 
conclusion that the evolution of the condensate wave function ip t should be described by the nonlinear 
Hartree equation (|2.9p with V replaced by bo 5, that is by the equation 

id t (ft = -Acpt + b (5 * \(pt\ 2 )pt = -Apt + bo\(pt\ 2 (ft ■ (3.16) 

Comparing with (|3.13p . we note that (|3.16p is characterized by a different coupling constant in front 
of the nonlinearity. The emergence of the scattering length in the Gross-Pitaevskii equation (|3.13p is 
the consequence of a subtle interplay between the iV-dependent interaction potential and the short 
scale correlation structure developed by the solution of the iV-particle Schrodinger equation ipN,t 
(as we will see, the correlation structure varies on lengths of the order 1/N, the same lengthscale 
characterizing the interaction potential). This remark implies that, in order to prove the convergence 
to the infinite hierarchy (where the coupling constant 87rao already appears), we will need to identify 
the singular correlation structure of ipN,t (which is then inherited by the marginal densities 7jy j)- 
This is one of the main difficulties in the proof of the convergence, which was completely absent 
in the analysis of mean-field systems presented in Section [2j we will discuss it in more details in 
Section 13.21 and in Section 13.31 

The presence of the correlation structure in ipN,t also affects the proof of a-priori bounds 

Tt(l-A 1 )...(l-A fc ) 7 £!t<C* (3.17) 

for all k > 1, t € H., for the limit points T^ = {7^ t }fc>i of the marginal densities Tjv,* = {l^t\k=i 
associated with 4>N,t ■ As in the case of a Coulomb potential discussed in Section 12.41 these a-priori 
bounds play a fundamental role because they allow us to restrict the proof of the uniqueness to a 
smaller class of densities. In Section \2A\ we derived the a-priori bounds for 7^ making use of the 
estimates (I2.46P which hold uniformly in N, for N sufficiently large (in turns, the bounds (|2.46p were 
obtained through the energy estimates of Proposition 12. 7p . In the present setting, however, bounds 
of the form 

Tr(l-A 1 )...(l-A k )^ t <C k (3.18) 

cannot hold uniformly in N, because of the short scale correlation structure developed by the solution 
of the Schrodinger equation. Remember in fact that the short scale structure varies on the length 
scale 1/N; therefore, when we take derivatives of 7Jy 4 as in (|3.18j) . we cannot expect to obtain bounds 
uniform in N (unless we take only one derivative, because of energy conservation). Although the 

marginals "y Nt , for large but finite N, do not satisfy the strong estimates (I3.18|) . it turns out that 

one can still prove the a-priori bound (|3.17p on the limit point 7^ ' t . This is indeed possible because 
in the weak limit N — ► 00, the singular short scale correlation structure characterizing the marginal 
densities j Nt disappears, producing limit points 7^ which are much more regular than the densities 

7jy- r Because of the absence of estimates of the form (|3.18p for j Nt , the proof of the a-priori bounds 

(k) 

for the limit points 7^ t requires completely new ideas with respect to what has been discussed in 
Section 12.41 we briefly discuss the most important ones in Section 13.41 

Finally, the singularity of the interaction potential strongly affects the proof of the uniqueness 
of the solution to the infinite hierarchy. In Section 12.41 the main idea to prove the uniqueness of 
the infinite hierarchy was to expand the solution in a Duhamel series and to control all Coulomb 
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potentials appearing in the expansion through Laplacians acting on appropriate variables and at the 
end to control the expectation of the Laplacians through the a-priori bounds (J2.46P on the densities 

(k) 

Too £. In this argument, it was very important that the Coulomb potential can be controlled by the 
kinetic energy, in the sense of the operator inequality 

-!-<C(l-A). (3.19) 

\x\ 

In the present setting, the Coulomb potential has to be replaced by a 5-function. In three dimensions, 
the (^-potential cannot be controlled by the kinetic energy. In other words, the bound 

S(x) <C(1-A) a 

is not true for a = 1; it only holds if a > 3/2 (in three dimensions, the L°° norm of a function 
can be controlled by the i7 Q -norm, only if a > 3/2). This observation implies that the a-priori 
bounds (|3.17p are not sufficient to conclude the proof of the uniqueness of the infinite hierarchy with 
delta-interaction (while similar bounds were enough to prove the uniqueness of the infinite hierarchy 
with Coulomb potential). Since it does not seem possible to improve the a-priori bounds to gain 
control of higher derivatives (one would need more than 3/2 derivatives per particle), we need new 
techniques to prove the uniqueness of the infinite hierarchy. We will briefly discuss these new methods 
in Section 13.51 

3.2 Convergence to the Infinite Hierarchy 

The goal of this section is to discuss the main ideas used to prove the next proposition which identifies 
limit points of the sequence Tjv,t = {In t\k=i as solutions to a certain infinite hierarchy of equations 
(this proposition replaces Proposition 12.31 which was stated for mean-field systems with bounded 
interaction potential). 

Proposition 3.2. Suppose that V > 0, with V(x) < C{x)~ a , for some a > 5, and for all x € M 3 . 
Assume that the sequence i/jn satisfies i3.11\) and the additional assumption \3. 1$ . Fix T > and 

let Yooj = {7oo,*}fc>i G 0fc>iC([°' r ]'^fc) be a limit point ofT N>t = {TN,t\k=i ( with respect to the 
product topology T pro( i defined in Section \2.3\) . Then Tqc^ is a solution to the infinite hierarchy 

k i 

TlV^hSo-SwE / dsU^{t- S )Tr k+l \6( Xj -x k+1 ),^+V] (3.20) 

with initial data <y^ Q = |^)( (/9 |® fc (see [2A3\) for the definition ofU^>). 

The detailed proof of this proposition can be found in |15|, Theorem 8.1] (for small interaction 
potential, see also [13], Theorem 7.1]). 

(k) 

To prove the proposition, we start by studying the time-evolution of the marginal densities *y Nt , 
which is governed by the BBGKY hierarchy. In integral form, the BBGKY hierarchy is given by 

k 



^\=U^\t) 1 %- l J2 I dsU^it-s) \v N ( Xi - Xj ),^ 
i<j Jo L 



„( fc ) 



(3.21) 



(N-k)f2 f dsU^(t-s)Tr k+1 \v N ( Xj -x k+1 ),^ 
~1 Jo L 
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(k+i) 



Assuming (by passing to an appropriate subsequence) that T^t — ► r^^ as N — > oo with respect to 
the product topology r pro( i introduced in Section [2T3l it is simple to prove that the l.h.s. and the first 
term on the r.h.s. of (|3.21|) converge, as N — ► oo, to the l.h.s. and, respectively, to the first term on 
the r.h.s. of (|3.20p . The second term on the r.h.s. of (|3.2ip . on the other hand, can be proven to 
vanish in the limit N —* oo (at least formally, this follows by the observation that the second term 
is smaller by a factor of N w.r.t. the third term). The fact that the second term on the r.h.s. of 
(j3.2ip is negligible in the limit N — > oo (compared with the third term) corresponds to the physical 
intuition that the interactions among the first k particles affect their time-evolution less than their 
interaction with the other (N — k) particles. 

To conclude the proof of Proposition 13.21 we only need to show that the third term on the r.h.s. 
of (|3.2ip converges, as N — ► oo, to the last term on the r.h.s. of (|3.20p . As already remarked in 
Section [3.14 this convergence relies critically on the correlation structure characterizing the (k + 1)- 
particle density -y Nt . A naive approach, based on the observation that (N — k)VN(xj — x k +i) ~ 
N 3 V(N(xj — Xfe+i)) ~ bo5(xj — Xfe+i) for large N, fails to explain the coupling constant in front of the 
last term on the r.h.s. of (|3.20p . The emergence of the scattering length can only be understood by 

(k-\-l) 

taking into account the correlation structure of ^y N t . Assuming for a moment that the correlations 
can be described, in good approximation, by the solution /jv to the zero-energy scattering equation 
(|3.5p . we can expect that, for large N, 



tnj ( x fc+i;4+i) - In(xj - xjH-ihSj^xfc+isx'fc+i) ( 3 - 22 ) 

in the region where Xj — x k+ \ is of the order 1/N (and all other variables are at larger distances). 
Assuming some regularity of the limit point 7^ t , and using (|3.4p , the approximation (13.221) imme- 
diately leads to 



Tr fc+ i(iV - k)V N (xj - x fc+ i)7^ 1) j(x fe ;x^) 

- / dx k+1 N 3 V(N(xj - x k+1 ))f(N(xj -x fc+ i))7^ 1) (x fc ,j; fc+ i;x / fc ,x/ c+ i) 

= J dy V{y)f{y)^ t l) (x k , Xj + -^; x' k , Xj + -|j (323) 

~ (J 'dyV(y)f(y)) j^^xj^Xj) 

= 87ra / dx k+ iS(xj - Xfc+i^t ( x fc>^fe+i; x fc>^fc+l) 



and thus explains the emergence of the scattering length on the r.h.s. of (|3.20p (note that the third 
term on r.h.s. of (|3.2ip is a commutator and thus produces two summands; in (13.230 we only consider 
one of these terms, the other can be handled analogously). This heuristic argument shows that in 
order to prove Proposition 13.21 we need to identify the short scale structure of the marginal densities 
and prove that it can be described by the function fjy as in (|3.22p . To this end we are going to 
use energy estimates. In [13] and [15j . we developed two different approaches to this problem. The 
first approach is simpler, but it only works for sufficiently small interaction potentials. The second 
approach is a little bit more involved, but it can be used for all potentials satisfying the assumptions 
of Theorem 13. 11 In the following we will focus on the first, simpler, approach; in the next subsection, 
we present the main ideas of the second approach. 

To measure the strength of the interaction potential V, we define the dimensionless constant 

p = sup \x\ 2 V{x) + I prV{x). (3.24) 

xem. 3 J \ x \ 
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Proposition 3.3. Assume that the potential V satisfies the conditions of Theorem \3.1[ and suppose 
that p > is sufficiently small. Then there exists C > such that 

2 



(^H 2 N rP)>CN 2 



dx 



ViVj 



^(x) 



f N (Xi -Xj) 



(3.25) 



for all i ^ j and for all ifi £ L\ 



,dx). 



This energy estimate, combined with the assumption (|3.14p on the initial wave function vp^, leads 
to the following a-priori bounds on the solution ifiwt = e~ lHNt ip]\f of the Schrodinger equation (|3.9p . 

Corollary 3.4. Assume that V satisfies the conditions of Theorem \3.1\ and suppose that p > is 
sufficiently small. Suppose that ipN satisfies \3. 1 0\) and ft3.14\ ). Then we have 

2 



I 



dx 



^N,t( x ) 



^Vj— 

jN(Xi ~ Xj) 



<C 



(3.26) 



(k) 

for all i ^ j, uniformly in N £ N and in t E R. Therefore, ifjpji denotes the k-particle marginal 
associated with ipN,t, we have, for every 1 <i,j < k with i ^ j, 



Tr(l-Ai){l-A 



(k) 



f N {Xi - Xj) JV ' t f N (x t - Xj) 



<C 



uniformly in N G N and in t £ 



Proof. Using (|3.25|) . the conservation of the energy along the time evolution, and the assumption 
(|3.14p on the initial wave function ^jy, we find 



dx 



ViV ._?Mx)_ 

3 fN{xi - Xj) 



< CN- 2 {i, N}t ,H^ N!t ) = CN- 2 (^ N ,H 2 N ^ N ) < C. 



U 



Remark that the a-priori bounds (J3.26P cannot hold true if we do not divide the solution i/jjsrt of 



the Schrodinger equation by /at(xj 
to check that 



). In fact, using that /jv(a?) — 1 — ao/(N\x\ + 1), it is simple 



fdx\V 2 f N (x) 



N. 



This implies that, if we replace V , iv,t( x )//iv( a; i 



by V"at(x) the integral in (|3.26p would be of order 



N. Only after removing the singular factor fN{x% —Xj) from ^7v,t( x ) we can obtain useful bounds on 
the regular part of the wave function (regular in the variable (xj — Xj)). These a-priori bounds allow 
us to identify the correlation structure of the wave function ipN,t and to show that, when Xi and 
Xj are close to each other, i^N,t( x ) can be approximated by the time independent correlation factor 



fN(Xi 



u 3Ji 



which varies on the length scale 1/iV, multiplied with a regular part (which only varies 



on scales of order one). In other words, the bounds (J3.26P establish a strong separation of scales for 
the solution ipN,t of the ./V-particle Schrodinger equation, and for its marginal densities; on length 
scales of order 1/N, ipN,t is characterized by a singular, time independent, short scale correlation 
structure described by the the solution /jy to the zero-energy scattering equation. On scales of order 
one, on the other hand, the wave function ipN,t is regular, and, as it follows from Theorem 13 .1\ it can 
be approximated, in an appropriate sense, by products of the solution to the time-dependent Gross- 
Pitaevskii equation. Remark that although the short-scale correlation structure is time independent, 
it still affects, in a non-trivial way, the time-evolution on length scales of order one (because it 
produces the scattering length in the Gross-Pitaevskii equation). 
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Proof of Proposition \3.3l We decompose the Hamiltonian (|3.8p as 

N 1 

Hn = 2. hj with hj = — Aj + — 2_^ Vn( 



2 

For an arbitrary permutation symmetric wave function i\) and for any fixed i ^ j, we have 

(V, H 2 N i>) = N(i/>, h 2 t iP) + N(N - l)(ij,, hihrf) > N(N - l)(ij,, hihjiP) . 
Using the positivity of the potential, we find 

(V, Hfap) > N(N - 1) U, f-Ai + ~V N (xi - Xj )\ (-Aj + ^V N ( Xi - ay) J A . (3.27) 

Next, we define </>(x) by ?/>(x) = /iv(^i — ay) ^>(x) (0 is well defined because /at (x) > for all x £ 
note that the definition of the function <p depends on the choice of i, j. Then 

1 \.r*..r„. .uu\_Aju , (A/jv)(xi - Xj) . . V/at(xj - ay), 



-A* (/jv(xi - xj)0(x)) = Ai^(x) + v ; J 7 V r^0(x) + -^ f V i( /.(x) . 



fN(Xi~Xj) fN(Xi-Xj) " f N (Xi-Xj) 

From f)3.3f) it follows that 

— — - — - ( -A, + -V N (xi - Xj) ) f N (xi - Xj)<j){x) = Li(f>(x) 

JN\ x i Xj) \ 2 

and analogously 



A j + T, V N{Xi ~ Xj) 1 f N (Xi - Xj)(f>(x) = Lj(f>(x) 



f N (xi -Xj) V 2 

where we defined 

T A , V^/Ar(xi-Xj) . . 

L £ = -A^ + 2— — J-S7e, for l = i,j- 

fN(Xi - Xj) 

Remark that, for I = i,j, the operator Li satisfies 

dxfx(xi-Xj) L £ (f)(x.) ip(x) = / dx fx(xi-Xj) 0(x) L £ V( X ) = / dx/^(x;-x-,) V £ 0(x) V^(x) . 
Therefore, from (|3.27p . we obtain 

(iP,H 2 N iP) > N(N - 1) Jdx f 2 N { Xi - Xj) L^(x) Lj 4>(x) 

= N(N - 1) J dx /^(x, - ay) V^(x) ViL, 0(x) 

= N(N - 1) / dx /^(xi - x,) V^(x) Lj V,0(x) 

. _ (3.28) 

+ N(N - 1) / dx /^(xi - xj) V^(x) [V^^x) 



JV(iV - 1) j dx /^(x 4 - xj) |Vj-Vi0(x) 



+ .V(.Y-J.) /dx/^fc-ay) (vi V , fN f iXi ^ ) V,,Xx)V,o(x). 

V JN{Xi-Xj) 
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To control the second term on the right hand side of the last equation we use bounds on the function 
/at, which can be derived from the zero energy scattering equation f)3.3() : 



l-Cp<f N (x)<l, \Vf N (x)\<Cf-, \V 2 f N (x)\<C-^ (3.29) 



, 19 

\x\ \x\ z 



for constants C independent of iV and of the potential V (recall the definition of the dimensionless 
constant p from (|3.24p ). Therefore, for p < 1, 



dxfaxi-Xj) V t v 7 T V*)Vi(x 



f N (Xi - Xj) 

< Cp I dx 1 —— |V^(x)| |V^(x) 



LX^ 7 I 



< Cp /dx -j (|V^(x)| 2 + |V^(x)| 2 ) 

J | X{ Xj I 

<Cp /dx|V i V j 0(x)| 2 



where we used Hardy inequality. Thus, from (|3.28p . and using again the first bound in (|3.29|) . we 
obtain 

(t/j,H%^) > N(N-l)(l-Cp) /dx|ViV^(x)| 2 

which implies (|3.25p . D 

Equipped with the a-priori bounds of Corollary 13.41 we can now come back to the problem of 
proving the convergence of the last term on the r.h.s. of (|3.21|) to the last term on the r.h.s. of 
(|2.6p . For simplicity, we consider the case k = 1, and we only discuss the term with the interaction 
potential on the left of the density (the commutator also has a term with the interaction on the right 
of the density, which can be handled analogously). After multiplying with a smooth one-particle 
observable j' 1 ' (a compact operator on L 2 (M 3 ), with sufficiently smooth kernel), we need to prove 
that 

Tr (U^{s - t)J®>) (n z V{N{ Xi - x 2 )h5S " 8W(xi - ^b^t) ^ 

as N —>■ oo. To this end we decompose the difference in several terms. We use the notation J t = 
U^'(t) J^ 1 ' , and, for a bounded function h{x) > with f dx h(x) = 1, we define h a {x) = a~ 3 h(a^ 1 x) 
for all a > 0. Then we have 

Tr (w«(s - t) jW) (n 3 V(N( Xi - x 2 ))7$ - 87ra 5( Xl - x 2 ) 7 g t ) 

= Tr J« N*V(N( Xl - X2 ))f(N( Xl - ^)) /( ^ (x ;_ X2)) 7% /(iV(xi 1 _ X2)) (/(iV(x 1 - *,)) - 1) 

+ Tr jffi (N*V(N( X1 - X2 ))f { N {xi - x 2 )) - S,a S( Xl - x 2 )) —-i-^ % —^^ 

+ 8na Tr J« (S( X1 - x 2 ) - h a { Xl - x 2 )) f{N{x \_ X2)) l { N\ f(N(x \_ X2)) 

+ 8vra Tr J« h a ( Xl - x 2 ) ( f^.^f ^ f^.^ " 7$) 

+ 8vra Tr J s a) t /i a (xi - x 2 ) (j^t ~ T^t, 



+ 87ra Tr Jf\ {h a {x\ - x 2 ) - S(xi - x 2 )) j£ t . 



(3.30) 
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The idea here is that in order to compare the ^-dependent potential N^V(N(xi — X2)) with the 

limiting <5-potential, we have to test it against a regular density (using an appropriate Poincare 

(2) 
inequality). For this reason, we first regularize the density *y Nt in the variable (x± — X2) dividing it 

by the correlation function fm(xi — x-z) on the left and the right (first term on the r.h.s. of the last 
equation). Using the regularity of f N ^{x\ — £2)7jVt/iV ( Xl ~~ X2 ) fr° m Corollary 13.44 we can then 
compare, in the regime of large N, the interaction potential with the delta-function (second term on 
the r.h.s.). At this point we are still not done, because, in order to remove the regularizing factors 
/jv ( x i ~~ x 2) (fourth term on the r.h.s. of (|3.3U|) ) and in order to replace the density ~/ Nt by its 

limit point 7^ (fifth term on the r.h.s. of (|3.30j) ). we need to test the density against a compact 

observable. For this reason, in the third term on the r.h.s. of (|3.30p . we replace the <5-function (which 

is of course not bounded) by the function h a which approximate the delta-function on the length 

scale a; it is important here that a is now decoupled from N. In the last term, after removing all 

(2) 
the N dependence, we go back to the 5-potential using the regularity of the limiting density 7^ t . 

To control the first and fourth term on the r.h.s. of (J3.30J) . we use the fact that 1 — /jv(xi — X2) — 
l/(N\x\ — X2I + 1) varies on a length scale of order 1/N. It follows that the first term converges 
to zero as N — > 00, as well as the fourth term, for every fixed a > 0. To estimate the second, the 
third and the last term, we make use of appropriate Poincare inequalities, combined with the result 
of Corollary 13.41 and, for the last term, of Proposition 13.71 (we present an example of a Poincare 
inequality, which can be used to estimate these terms in Appendix |A|) . It follows that the second 
term converges to zero as N — > 00, and that the third and the fifth terms converge to zero as a — > 0, 

uniformly in N. Finally, the fifth term on the r.h.s. of (13.301) converges to zero as N — ► 00, for every 

(2) (2) 

fixed a; this follows from the assumption that ~/ Nt — ► 7^ t as N — ► 00 with respect to the weak* 

topology (some additional work has to be done here, because the operator J s _ t h a {x\ — X2) is not 
compact). Therefore, if we first fix a > and let N — > 00 and then we let a — > all terms on the 
r.h.s. of (|3.30p converge to zero; this concludes the proof of Proposition 13.21 

3.3 Convergence for Large Interaction Potentials 

As pointed out in Section 13.21 the energy estimate given in Proposition 13.31 which was a crucial 
ingredient for the proof of Proposition 13. 2L only holds for sufficiently small potentials (for sufficiently 
small values of the parameter p defined in (|3.24p ). For large potentials, we need a different approach. 
The new technique, developed in [15] , is based on the use of the wave operator associated with the 
one-particle Hamiltonian fjjv = —A + (1/2)Vn, defined through the strong limit 

W N = s- lim e ii)Nt e iAt . (3.31) 

Under the assumptions of Theorem 13. II on the potential V, it is simple to show that the limit (|3.31[) 
exists, that the wave operator Wn is complete, in the sense that 

WZ l = W* N = s - lim e-^e-^Nt f 

t — >oo 

and that it satisfies the intertwining relation 

W^[)W N = -A. (3.32) 

It is also important to observe that the wave operator Wn is related by simple scaling to the wave 
operator W associated with the one-particle Hamiltonian f) = — A + (1/2) V (and defined analogously 
to (|3.3ip ). In fact, if Wn(x;x') and W(x;x') denote the kernels of Wn and, respectively, of W, we 
have 

W N (x;x') = N 3 W(Nx;Nx') and W N (x;x') = N 3 W*(Nx;Nx') . 

31 



In particular this implies that the norm of Wn, as an operator from L P (M. ) to L P (M ), for arbitrary 
1 < p < oo, is independent of N. From the work of Yajima, see |321 133] . we know that, under the 
conditions on V assumed in Theorem 13. 1\ W is a bounded operator from L P (M 3 ) to L P (M 3 ), for all 
1 < p < oo. Therefore 

|| Wn\\lp->lp = \\W\\lp^lp < oo for all 1 < p < oo . 

In the following we will denote by Wn u j\ the wave operator Wn acting only on the relative variable 
Xj — Xi. In other words, the action of Wn u j) on a iV-particle wave function tpN G L 2 (R 3N ) is given 

by 

// T* ' I T* ■ 1) T™ ■ I IT ~ 11 \ 

dv W N (xj-Xi;v)ipN ( xx,..., - - J +-,..., - - J - -, ...,xjv) (3-33) 

if j < i (the formula for i > j is similar). Similarly, we define W N u •-,. Using the wave operator we 
have the following energy estimate, which replaces Proposition 13.31 and whose proof can be found in 
[HI Proposition 5.2]. 

Proposition 3.5. Suppose V > 0, V G L 1 ^ 3 ) n L 2 (M 3 ) and V(x) = V(-x) for all x G M 3 . T/ien 
we have, for every i ^ j, 

Wn,H n tPn) > CN 2 J dx | (V< • V,) W^ i(i)j) ^7v| 2 • (3.34) 

From Proposition 13.51 we obtain immediately an a-priori bound on ipN,t and on its marginal 
densities. 

Corollary 3.6. Assume that V satisfies the conditions of Theorem \3.1\ Suppose that ipN satisfies 
\3.10\) and ^3.1$ - Then we have, for all i ^= j, 

/dx |(Vi • Vj) Wh (iJ) il>N,t(x)\ 2 ^ c ( 3 - 35 ) 

(k) 

uniformly in N G N and t G K. Therefore, if "fj^i denote the k-particle marginal associated with 
ipN,t, we have, for every 1 < i,j < k with i / j, 

Tr ((V, • V,) 2 - A, - A d + 1) W Njiid) ^ t WN, iid) < C 

uniformly in N G N and in t G R. 

The philosophy of the bounds (13,350 and (13.260 is the same; first we have to regularize the wave 
function ipN,tj and then we can prove useful bounds on its derivatives. There are however important 
differences. In (13.260 we regularized ipN,t in position space, by factoring out the short scale correlation 
structure /at(xj — Xj). In (|3.35p . instead, we regularize ipN,t applying the wave operator W N a-y 
Another important difference is that (|3.35p is weaker than (13.260 ; in fact, (13.350 only gives a control on 
the combination ^2 a= id Xia d x - a , while (|3.26p controls d Xia d x . g for all 1 < a, (3 < 3. The weakness 
of the bound (J3.35P makes the proof of the convergence more difficult. In particular we have to 
establish new Poincare inequalities, which only require control of the inner product Vj • Vj. It turns 
out that the weaker control provided by (J3.35P is still enough to conclude the proof of convergence 
to the infinite hierarchy (Proposition I3.2P , For more details, see [15j. Section 8]. 
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3.4 A-Priori Estimates on Limit Points T^j 

In this section we present some of the arguments involved in the proof of the a-priori bounds (|3.17p . 

Proposition 3.7. Assume that V satisfies the conditions of Theorem \3.1\ Suppose that i/jn satisfies 
i3.10\) and {3.1$ . Let T^j = {7oo t }fc>i 6 ©fc>i C([0, T},C\) be a limit point of the sequence 

^N,t = {TyvtlfcLi with respect to the product topology T pro d defined in Section \2. 3\ Then 7^ ' t > and 
there exists a constant C such that 

Tr(l-A 1 )...(l-A k )^] t <C k (3.36) 

for allk>l andt G [0,T\. 

The main difficulty in proving Proposition 13.71 is the fact that the estimate (|3.36p does not hold 

(k) (k) 

true if we replace 7^ with the marginal density j N \. More precisely, 

Tr (1 - Ax) ... (1 - A fc ) 7 $ < C k (3.37) 

cannot hold true with a constant C independent of N. In fact, for finite N and k > 1, the /c-particle 

(k) 

density ^ N t still contains the singular short scale correlation structure. For example, when particle 
one and particle two are very close to each other (at distances of order 1/N), we can expect the 
two-particle density to be approximately given by 

T N 't(^2,x 2 ) - const In(xi - x 2 )Jn(x 1 - x' 2 ) 

(the constant part takes into account factors which vary on larger scales). It is then simple to check 
that 

Tr(l-A 1 )(l-A 2 ) 7 g ) t ^iV. 

Only after taking the weak limit N — > 00, the short scale correlation structure disappears (because 
it varies on a length scale of order 1/N), and one can hope to prove bounds like (|3.36[) . 

To overcome this problem, we cutoff the wave function ipN,t when two or more particles come 
at distances smaller than some intermediate length scale £, with iV -1 <C £ <C 1 (more precisely, the 
cutoff will be effective only when one or more particles come close to one of the variable Xj over 
which we want to take derivatives). For fixed j = 1, . . . , N, we define 9j 6 C°°(lR 3Ar ) such that 

. , . J 1 if I x i — Xj I ^> £ for all i 7^ j 
3 I if there exists i 7^ j with \xi — Xj\ < £ 



It is important, for our analysis, that 9j controls its derivatives (in the sense that, for example, 

|Vj0j| < C£~ 1 6- ); for this reason we cannot use standard compactly supported cutoffs. Instead 
we have to construct appropriate functions which decay exponentially when particles come close 
together (the prototype of such function is 9{x) = exp(—y/(x/£) 2 + 1)). Making use of the functions 
9j(x), we prove the following higher order energy estimates. 

Proposition 3.8. Choose £ <C 1 such that N£ 2 3> 1. Suppose that a is small enough. Then there 
exist constants C\ and C2 such that, for any ip G L 2 (R. ), 

W>, (H N + CiAOV) > C 2 N k J dx 0i (x) . . . fc _i(x) |Vi . . . V fc V(x)| 2 . (3.38) 
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The meaning of the bound (|3.38[) is clear. The L 2 -norm of the fc-th derivative Vi . . . Vfc^ can 
be controlled by the expectation of the k-th power of the energy per particle, if we restrict the 
integration domain to regions where the first (k — 1) particles are "isolated" (in the sense that there 
is no particle at distances smaller than £ from xi,X2,---, Xk-l)- 

Note that we can allow one "free derivative"; in (J3.38P we take the derivative over x^ although 
there is no cutoff 0&(x). The reason is that the correlation structure becomes singular, in the L 2 - 
sense, only when we derive it twice (if one uses the zero energy solution /jy introduced in (|3,3p 
to describe the correlations, this can be seen by observing that |V/jv(x)| < l/\x\, which is locally 
square integrable). Remark that the condition N£ 2 3> 1 is necessary to control the error due to the 
localization of the kinetic energy on distances of order £. The proof of Proposition 13.81 is based on 
induction over k; for details see Section 7 in 1151 . 



From the estimates (I3.38D , using the preservation of the expectation of H^ along the time evo- 
lution and the condition (|3.14p . we obtain the following bounds for the solution i/jN,t = e~ lHNt ijJN of 
the Schrodinger equation (|3.9p . 

dx 0x(x) . . . fc _i(x) | Vi . . . V fc ^,t(x)| 2 < C k 



I 



uniformly in N and t, and for all k > 1. Translating these bounds in the language of the density 
matrix 7jv,t, we obtain 

Tr 6 1 . . . fc _iV! . . . V fc7iV ,tVl ...V* k <C k . (3.39) 

The idea now is to use the freedom in the choice of the cutoff length £. If we fix the position of all 
particles but Xj, it is clear that the cutoff 8j is effective at most in a volume of the order N£ 3 . If 
we choose £ such that N£ 3 — ► as N — ► oo (which is of course compatible with the condition that 
N£ 2 3> 1), we can expect that, in the limit of large N, the cutoff becomes negligible. This approach 
yields in fact the desired results; starting from (|3.39p . and choosing £ such that N£ 3 < 1, we can 
complete the proof of Proposition 13.71 (see Proposition 6.3 in [13] for more details). 

3.5 Uniqueness of the Solution to the Infinite Hierarchy 

To complete the proof of Theorem 13.11 we have to prove the uniqueness of the solution to the infinite 
hierarchy (]3.20p in the class of densities satisfying the a-priori bounds Q3.36p . Remark that the 
uniqueness of the infinite hierarchy (|3.20p . in a different class of densities, was recently proven by 
Klainerman and Machedon in [25]. The proof proposed by Klainerman and Machedon is simpler 
than the proof of Proposition 13.91 which we discuss below. Unfortunately, the result of [25] cannot be 
applied to the proof of Theorem l3.lt because it is not yet clear whether limit points of the sequence 
of marginal densities IV,* = {7jvt}fc^=l fit into the class of densities for which uniqueness is proven. 

Proposition 3.9. Fix T > and T = {j^}k>i & ©fc>i C\- Then there exists at most one solution 

T t = {^ k) } k >i G ©C([0,T],£ fc ) of the infinite hierarchy ETWj) with r t=0 = V, such that 7 f } > 
is symmetric with respect to permutations, and 

Tr(l-A 1 )...(l-A fc ) 7 ( fc ) <C k (3.40) 

for all k > 1 and all t G [0, T]. 



In this section we briefly explain some of the main steps involved in the proof of Proposition 13.9 
the details can be found in [12] [Section 9] . 
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To shorten the notation, we write the infinite hierarchy (|3.2Q[) in the form 

lt = U ik \t)-f + [ ds U^it - s) B^^f +1) , 



(3.41) 



where Zv ■ '(t) denotes the free evolution of k particles 



U {k \t)^ 



itH% 1 ^ 1 {k) e -itY. k j=l ^ 



and the collision operator B^ k > maps (k + l)-particle operators into A;-particle operators according to 

k 



i? (fe) 7 (fc+1) = -8ivra ^Tr fc+1 [s( Xj 

3=1 



Xk+l),l 



(fc+1) 



(3.42) 



The map B^ ' is defined in analogy to Section [2] in particular the kernel of B^ >^ +1 > is given by 
the expression on the r.h.s. of (|2.1ip . with V(x) replaced by 8irao5(x). 

Iterating (|3.4ip n times we obtain the Duhamel type series 



.,(*) 



n-1 



(fc) 



7f ) =^)(t) 7 r+x:ci+< 



:,t 



(3.43) 



m=l 



with 






t 

dsi . . 
o 

k fc+1 

EE 

31=1 32=1 j m =l 



ds m U^it - s 1 )B { - k) U { - k+1 \s 1 - s 2 )B 



(fc+i) 



B (k+ m -i) u {k+ m )^ Sm) ^ 



(fc+m) 



k+m „t 



ds m U {k \t - si)Tr fc+ i 5{x h - x k+ i) 



U^ k+l \ Sl - s 2 )Tr k+2 \s(x h - x k+2 ), . . . Tr fc+m \d(x jm - x k+m )M {k+m \s m )l { k+m) 



(3.44) 



and the error term 

Vnl = / ds i / ds 2 ... 

Jo Jo Jo 



ds n U^ k \t - Sl )B^U^ k+l \ Sl - S 2 )B( k+ V . . . £(fc+n-l) 7 (fc+m) _ 



(3-45) 

Note that the error term (|3.45p has exactly the same form as the terms in (|3.44j) . with the only 
difference that the last free evolution is replaced by the full evolution js n • 

To prove the uniqueness of the infinite hierarchy, it is enough to prove that the fully expanded 
terms (|3.44[) are well-defined and that the error term (J3.45P converges to zero as n — ► oo (in some 
norm, or even after testing it against a sufficiently large class of smooth observables) . The main 
problem here is that the delta function in the collision operator B' ' cannot be controlled by the 
kinetic energy (in the sense that, in three dimensions, the operator inequality 5{x) < C(l — A) does 
not hold true). For this reason, the a-priori estimates (|3.40p are not sufficient to show that (J3.45P 
converges to zero, as n — > oo. Instead, we have to make use of the smoothing effects of the free 
evolutions U^ k+3 '(sj — Sj+i) in (|3.45p (in a similar way, Stricharzt estimates are used to prove the 
well-posedness of nonlinear Schrodinger equations). To this end, we rewrite each term in the series 
(|3.43|) as a sum of contributions associated with certain Feynman graphs, and then we prove the 
convergence of the Duhamel expansion by controlling each contribution separately. 
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Figure 1: A Feynman graph in T m ,k and its two types of vertices 



The details of the diagrammatic expansion can be found in |X2(, Section 9]. Here we only sketch 
the main ideas. We start by considering the term Q n \ in (|3.44p . After multiplying it with a compact 



fc-particle observable J^ k > and taking the trace, we expand the result as 



Tr J«e (fc) 



m,t 



E ^ 



(3.46) 



Aer n 



where K\j is the contribution associated with the Feynman graph A. Here J- m ,k denotes the set of 
all graphs consisting of 2k disjoint, paired, oriented, and rooted trees with m vertices. An example 
of a graph in J- m ,k is drawn in Figure HJ Each vertex has one of the two forms drawn in Figure [U 
with one "father"-edge on the left (closer to the root of the tree) and three "son"-edges on the right. 
One of the son edge is marked (the one drawn on the same level as the father edge; the other two 
son edges are drawn below). Graphs in J- m k have 2k + 3m edges, 2k roots (the edges on the very 
left), and 2k + 2m leaves (the edges on the very right). It is possible to show that the number of 
different graphs in T m ^ is bounded by 2 4m+k . 

The particular form of the graphs in T m) k is due to the quantum mechanical nature of the 
expansion; the presence of a commutator in the collision operator (|3.42p implies that, for every 
^( fc +i) i n (|3,44p , we can choose whether to write the interaction on the left or on the right of the 
density. When we draw the corresponding vertex in a graph in J- m ,ki we have to choose whether to 
attach it on the incoming or on the outgoing edge. 

Graphs in J- m k are characterized by a natural partial ordering among the vertices (v -< v' if 
the vertex v is on the path from v' to the roots); there is, however, no total ordering. The absence 
of total ordering among the vertices is the consequence of a rearrangement of the summands on 
the r.h.s. of (|3.44j) ; by removing the order between times associated with non-ordered vertices we 
substantially reduce the number of terms in the expansion. In fact, while (|3.44j) contains {m + k)\/k\ 
summands, in (|3.46p we are only summing over at most 2 4m+fc contributions. The price we have to 
pay is that the apparent gain of a factor 1/m! due to the ordering of the time integrals in (|3.44|) is 
lost in the new expansion (|3.46p . However, since we want to use the time integrations to smooth 
out singularities it seems quite difficult to make use of this factor 1/m!. In fact, we find that the 
expansion (13.460 is better suited for analyzing the cumulative space-time smoothing effects of the 
multiple free evolutions than (|3,44p . 
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Because of the pairing of the 2k trees, there is a natural pairing between the 2k roots of the 
graph. Moreover, it is also possible to define a natural pairing of the leaves of the graph (this is 
evident in Figure [Q ; two leaves t\ and £2 are paired if there exists an edge e\ on the path from t\ 
back to the roots, and an edge ei on the path from £2 to the roots, such that e\ and e2 are the two 
unmarked son-edges of the same vertex (or, in case there is no unmarked sons in the path from t\ 
and £2 to the roots, if the two roots connected to t\ and £2 are paired). 

For A G ^F m ,kj we denote by E{A), V(A), R{A) and L{A) the set of all edges, vertices, roots 
and, respectively, leaves in the graph A. For every edge e £ E{A), we introduce a three-dimensional 
momentum variable p e and a one-dimensional frequency variable a e . Then, denoting by 7q and 

by j( k > the kernels of the density 7q m> and of the observable J^ k > in Fourier space, the contribution 
K\,t in (|3.46|) is given by 

Wn^ n*(s>.W2>.) 

J ee£(A) e ^ e e ^ e veV(A) \e£v / \e£v / 

X exp -it ^ T e( a e +ir e fM e ) J^ {{Pe}e£R(A)) 7o {{Pe}e(LL{A)) • 

\ eeR(A) / 

(3.47) 

Here r e = ±1, according to the orientation of the edge e. We observe from (|3,47|) that the momenta 
of the roots of A are the variables of the kernel of J' ' , while the momenta of the leaves of A are the 
variables of the kernel of 7q (this also explain why roots and leaves of A need to be paired). 

The denominators (a e — pl+iT e n e )~ l are called propagators; they correspond to the free evolutions 
in the expansion (|3.44j) and they enter the expression (|3.47[) through the formula 

poo it(a+i(i) 

e ttp = / da 



a — p 2 + ifx 



(here and in (|3.4T|> the measure da is defined by da = 6! a /{2m) where d'a is the Lebesgue measure 



on 



The regularization factors // e in (|3.47|) have to be chosen such that //father = ^2 e = son ^e a * ever y 
vertex. The delta-functions in (J3.47P express momentum and frequency conservation (the sum over 
e £ v denotes the sum over all edges adjacent to the vertex v; here ±a e = a e if the edge points 
towards the vertex, while ±a e = — a e if the edge points out of the vertex, and analogously for ±p e ). 

An analogous expansion can be obtained for the error term rj^l in (I3.45h . The problem now 
is to analyze the integral (|3.4T|) (and the corresponding integral for the error term). Through an 
appropriate choice of the regularization factors \i e one can extract the time dependence of K\ >t and 
show that 

\K A , t \<c k+m t^ f n 7^!% n ^e^We^I 

J e eE(T) { ° e Pe) v&V(T) \eev J \eev J (3.48) 

x I J^ ({Pe} e eR(r)) | \t +m) ({Pe} e6i( r)) | 

where we introduced the notation (x) = (1 + x 2 ) 1 ' 2 . 

Because of the singularity of the interaction at zero, we may be faced here with an ultraviolet 
problem; we have to show that all integrations in (|3.48[) are finite in the regime of large momenta and 
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Figure 2: Integration scheme: a typical vertex 

large frequency. Because of (|3.40p . we know that the kernel 7q m ({pe}eeL(A)) m (|3.48|) provides 
decay in the momenta of the leaves. From (|3.40p we have, in momentum space, 

dp 1 ... dp n (pi + 1) . . . (pi + l)^ n) (pi, . . . ,p n ;pi, ■ ■ ■ ,p n ) < C n 

for all n > 1. Heuristically, this suggests that 

\t +m \{Pe}eeL { A))\ < ft {Pe)^' 2 , (3.49) 

eeL(A) 

where (p) = (1 + p 2 ) 1 ' 2 . Using this decay in the momenta of the leaves and the decay of the 
propagators (a e — p 2 )~ l ,e £ -^(A), we can prove the finiteness of all the momentum and frequency 
integrals in (|3.47p . On the heuristic level, this can be seen using a simple power counting argument. 
Fix k>1, and cutoff all momenta \p e \ > n and all frequencies \a e \ > k 2 . Each p e -integral scales 
then as k 3 , and each a e -integral scales as k 2 . Since we have 2k + 3m edges in A, we have 2k + 
3m momentum- and frequency integrations. However, because of the m delta functions (due to 
momentum and frequency conservation), we effectively only have to perform 2k + 2m momentum- 
and frequency-integrations. Therefore the whole integral in (|3.47p carries a volume factor of the 
order K 5 ( 2fc + 2m ) = K iofc+iOm_ ]\j OWj s i nce there are 2k + 2m leaves in the graph A, the estimate (I3.49P 
guarantees a decay of the order K - 5 / 2 ( 2fc + 2m ) = K -5k-5m^ r-j-^g 2k + 3m propagators, on the other 
hand, provide a decay of the order K ~ 2 ( 2fc + 3m ) = K - 4fc_6m . Choosing the observable J^ so that jW 
decays sufficiently fast at infinity, we can also gain an additional decay K~ ek . Since 

^lOk+lOm — 5k— 5m— 4k— 6m- 6k — m— 5k ^, -i 

rv ' r\i — rv <$^ 1 

for k > 1, we can expect (|3.47p to converge in the large momentum and large frequency regime. 
Remark the importance of the decay provided by the free evolution (through the propagators); 
without making use of it, we would not be able to prove the uniqueness of the infinite hierarchy. 

This heuristic argument is clearly far from rigorous. To obtain a rigorous proof, we use an 
integration scheme dictated by the structure of the graph A; we start by integrating the momenta 
and the frequency of the leaves (for which (J3.49J) provides sufficient decay). The point here is that 
when we perform the integrations over the momenta of the leaves we have to propagate the decay 
to the next edges on the left. We move iteratively from the right to the left of the graph, until we 
reach the roots; at every step we integrate the frequencies and momenta of the son edges of a fixed 
vertex and as a result we obtain decay in the momentum of the father edge. When we reach the 
roots, we use the decay of the kernel j' > to complete the integration scheme. In a typical step, we 
consider a vertex as the one drawn in Figure [2] and we assume to have decay in the momenta of the 
three son-edges, in the form \p e \~ , e = u,d,w (for some 2 < A < 5/2). Then we integrate over 
the frequencies a u ,ad,a w and the momenta p u -,PdiPw of the son-edges and as a result we obtain a 
decaying factor |p r | in the momentum of the father edge. In other words, we prove bounds of the 
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form 



/ 



da u da d da w dp u dp d dp w 5(a r = a u + a d - a w )5(p r = p u + p d - p w ) const 
\Pu\ x \Pd\ x \p w \ x (au-pl}(a d -p d }(a w -pl) \p r \ x ' 



Power counting implies that (|3.50p can only be correct if A > 2. On the other hand, to start the 
integration scheme we need A < 5/2 (from (J3.49P this is the decay in the momenta of the leaves, 
obtained from the a-priori estimates). It turns out that, choosing A = 2 + e for a sufficiently 
small e > 0, (|3.50p can be made precise, and the integration scheme can be completed. 

After integrating all the frequency and momentum variables, from (J3.48P we obtain that 

\K A ,t\ <C fc+m t m/4 

for every A E F m ,k- Since the number of diagrams in T m ,k is bounded by C k+m , it follows immediately 
that 

Tr j( fc ) ciSt < C k+m t m / 4 . 

Note that, from (|3.44j) . one may expect £^/ f to be proportional to t m . The reason why we only get 
a bound proportional to t" 1 ' 4 is that we effectively use part of the time integration to control the 
singularity of the potentials. 

The only property of 7q used in the analysis of (J3.47H is the estimate (J3.40J) . which provides 

the necessary decay in the momenta of the leaves. Since the a-priori bound (|3.40p hold uniformly 

(k) 

in time, we can use a similar argument to bound the contribution arising from the error term r] n [ 

in (J3.45P (as explained above, also ry nt can be expanded analogously to (|3.46p . with contributions 
associated to Feynman graphs similar to (J3.47P : the difference, of course, is that these contributions 
will depend on 7^ for all s G [0,t], while (|3.47p only depends on the initial data). We get 

TrJ^rQ. <C fc+n t n/4 . (3.51) 

This bound immediately implies the uniqueness. In fact, given two solutions Fit = {l\ t }fe>i an d 

^2,t = {l2t }k>i °f the infinite hierarchy (|3.4ip . both satisfying the a-priori bounds (|3.40p and with 
the same initial data, we can expand both in a Duhamel series of order n as in (J3.43J) . If we fix 
k > 1, and consider the difference between 7} t and 7^, all terms (|3.44p cancel out because they 
only depend on the initial data. Therefore, from (I3.5ip . we immediately obtain that, for arbitrary 
(sufficiently smooth) compact A;-particle operators JW, 



Tr J^ ( 7 g> - 7 g } ) I < 2 C k+n t n ' A 



Since it is independent of n, the left side has to vanish for all t < 1/(2C) 4 . This proves uniqueness 
for short times. But then, since the a-priori bounds hold uniformly in time, the argument can be 
repeated to prove uniqueness for all times. 

3.6 Other Microscopic Models Leading to the Nonlinear Schrodinger Equation 

As discussed in Section [3TTI the strategy used to prove Theorem l3.1l is dictated by the formal similarity 
with the mean-field systems discussed in Section [2j from (I3.15P the Hamiltonian characterizing dilute 
Bose gases in the Gross-Pitaevskii scaling can be formally interpreted as a mean field Hamiltonian 
with an A^-dependent potential converging to a delta-function as N — > 00 (the physics described by 
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the two models is however completely different). The choice of the iV-dependent potential Vn(x) = 
N 2 V(Nx) in the Gross-Pitaevskii scaling is, of course, not the only choice for which the formal 
identification with a mean-field model is possible. For arbitrary (3 > 0, we can for example define 
the iV-particle Hamiltonian 



Hsji = £ -A; + 1 £ N^V{N^Xi - x 3 )) 



N 

j = l Kj 

acting on the Hilbert space L^(R ). The Hamiltonian (|3.8|) is recovered by choosing (3 = 1. 
For < (3 < 1, the potential iV 3/3 y(iV^2;) still converges to a delta-function as iV — > oo, but 
the convergence is slower. This fact has important consequences for the macroscopic dynamics; it 
turns out, in fact, that for < (3 < 1 the correlation structure developed by the evolved wave 
function ipN,p,t = e~ %HN ^ 1 varies on much shorter length scales compared with the length scale N~P 
characterizing the potential. Therefore, for < (3 < 1, the time evolution of the condensate wave 
function is still governed by a cubic nonlinear Schrodinger equation; this time, however, the coupling 
constant in front of the nonlinearity is given by bo = J V instead of 87rao (recall that the emergence 
of the scattering length in the Gross-Pitaevskii equation was a consequence of the interplay between 
the correlation structure in the many body wave function and the interaction potential; since, for 
< (3 < 1, the potential and the correlation structure vary on different length scales, this interplay 
is suppressed). The following theorem can be proven using the techniques developed in [15] (the 
statement for < (3 < 1/2 was proven in [12] ; in [13] . the whole range < (3 < 1 was covered, but 
only for sufficiently small potentials). 



Theorem 3.10. Suppose that V > satisfies the same assumption as in Theorem \3.1\ and assume 
that < f3 < 1. Let ifj N (x) = Hf =1 <p(.Xj), for some cp G H 1 ^ 3 ) and ^ N>t = e~ iH ^ Nt Tp N with the 
mean-field Hamiltonian 

N N 

H P , N = 5>A,- + ±Y,N 3Pv (N ( '(x i -x j )). 

j=l i<j 

Then, for every fixed k > 1 and t 6l, we have 

7$ -MM* 

as N — ► oo, where (ft is the solution to the nonlinear Schrodinger equation 

id t <pt = -Av? t + o-\Lp t \ 2 if t 
with initial data (ft=o = P and with 

8irao if (3 = 1 



" 'b i/0</5< 1 

4 Rate of Convergence towards Mean-Field Dynamics 

From the results of Section [23 we obtain that, for every fixed t S R, and for every fixed k E N, 

(k) 

where ^ Nt is the A;-particle density associated with the solution ipN,t of the iV-particle Schrodinger 
equation, and (ft is the solution of the Hartree equation. From Theorem 12. II and Theorem 12.51 we do 
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not get any information about the rate of convergence of j Nt to \ift)((pt\ ■ We only know that the 

difference j Nt — \(ft)(<pt\® k converges to zero, but we do not know how fast. Also in Section we 
do not obtain any information about the rate of convergence of the iV-particle Schrodinger evolution 
towards the Gross-Pitaevskii equation. This is not only a question of academic interest; in order to 
apply these results to physically relevant situations, bounds on the error are essential. 

For bounded potential, (J2.25P implies (specializing to the case k = 1) that 



Tr 



(D 



1% t -\<Pt)(<Pt\ <CN~ X (4.1) 



for sufficiently short times < t < to = V(8||V||). Iterating the argument leading to (|4.ip to obtain 
estimates valid for larger times, it is possible to derive bounds of the form 



Tr 



.X 1 ) 



7$-|pt>tel|<^=r- (4-2) 

Although (|4.2p shows that, for every fixed t > 0, the difference j Nt — \(pt)(ipt\ converges to zero, it is 
not very useful in applications because it deteriorates too fast; one would like to find bounds on the 
difference ~/ Nt — \(pt)((pt\ which are of the same order in iV for every fixed time. 

In [28], a joint work with I. Rodnianski, we obtain such bounds for mean-field systems with 
potential having at most a Coulomb singularity; the problem of obtaining error estimates for the 
Gross-Pitaevskii dynamics is still open. To prove such bounds, we do not make use of the BBGKY 
hierarchy. Instead, we use an approach, introduced by Hepp in [21] and extended by Ginibre and 
Velo in [20], based on embedding the iV-body Schrodinger system into the second quantized Fock- 
space representation and on the use of coherent states as initial data (coherent states do not have a 
fixed number of particles; this is what makes the use of a Fock-space representation necessary). The 
Hartree dynamics emerges as the main component of the evolution of coherent states in the mean 
field limit. To obtain bounds on the difference j Nt — \<pt)(<Pt\ f° r initial data describing coherent 
states, it is therefore enough to study the fluctuations around the Hartree dynamics, and to prove 
that, in an appropriate sense, they are small. Since factorized iV-particle wave functions can be 
written as appropriate linear combinations of coherent states, the estimates for coherent initial data 
can be translated into bounds for factorized initial data. Using these techniques, we prove in [28] 
the following theorem. We focus in this section on three dimensional systems, which are the most 
interesting from the point of view of physics; most of the results however can be extended to dimension 

Theorem 4.1. Suppose that there exists D > such that the operator inequality 

V 2 (x) < D (1 - A*) (4.3) 

holds true. Let 

N 

^(x) = n^j)' 

for some (p £ H l (M. 3 ) with \\ip\\ = 1. Denote by ipN,t = e~ tHNt i/j]\[ the solution to the Schrodinger 
equation fj7j|) with initial data ipN,o = ^N , and let j Nt be the one-particle density associated with 
ipN,t- Then there exist constants C,K, depending only on the H 1 norm of (f and on the constant D 
on the r.h.s. of ( f^.3[ ) such that 

C 



Tr 
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for every tGl and every JVgN. Here (ft is the solution to the nonlinear Hartree equation 

idtipt = -Atp t + (V * \(p t \ 2 )vt (4.5) 

with initial data (ft=o = <*P- 



Remarks. Condition (|4.3p is in particular satisfied by bounded potentials and by potentials 
with an attractive or repulsive Coulomb singularity. Theorem 14.11 implies therefore Theorem 12.11 and 
Theorem 12.51 Note, moreover, that the decay of the order N~ 1 ' 2 on the r.h.s. of (|4.4|) is not expected 
to be optimal. In fact, for initially coherent states we obtain in Theorem 14.41 the expected decay of 
the order \/N for every fixed time t E R; unfortunately, when factorized initial data are expressed as 
a superposition of coherent states, part of the decay is lost (note, however, that for a certain class of 
bounded potential a decay of the order iV -1 for factorized initial data has recently been established 
in HDD. 



4.1 Fock-Space Representation 

We define the bosonic Fock space over L 2 (M. s ,dx) as the Hilbert space 

T = L 2 (M 3 , &x)®° n = C L 2 s (M. 3n , dxi . . . dx r 



n>0 n>l 



where we put L 2 (M 3 )® s0 = C. Vectors in T are sequences i/j = {'0 }n>o of re-particle wave functions 
ipw £ L 2 (M 3 ™). The scalar product on T is defined by 

(iplifo) = 5^m '^2 )X3(R3«) = Vl -02 +5Z / dx l ■ ■ ■ dx n1pl (Xl, ■ ■ ■ , Xn)^' (X!, . . . , X n ) . 
n>0 n>l J 

An N particle state with wave function tpjy is described on T by the sequence {ip^ n '}n>o where 
ip( n ) = o for all re ^ N and t/' ' = ipN- The vector {1, 0, 0, ... } £ T is called the vacuum, and will 
be denoted by Q. 

On J-, we define the number of particles operator J\f, by (Afipp n > = nip^ n '. Eigenvectors of J\f 
are vectors of the form {0, . . . , 0, •i/A" 1 ) , 0, . . . } with a fixed number of particles m. For /ei 2 ' 
we also define the creation operator a*(f) and the annihilation operator a(f) on T by 

1 n 
(a*(/)^) (n) (xi,...,x n ) = —=y^ f{xj)^ n ~ l \xi, ..., Xj-i,Xj+x, ...,x n ) 

3=1 

(a(/)^) (n) (xi, ... ,x n ) = V^+T [dx J{x)^ n+l Hx,x 1 , ...,x n ). 



The operators a*(f) and a(f) are unbounded, densely defined and closed; a*(f) creates a particle 
with wave function /, a(f) annihilates it. It is simple to check that, for arbitrary n > 1, 

(°W n = {o,... ,o, /•»,(>,. ..}. 



The creation operator o*(/) is the adjoint of the annihilation operator a(f) (note that by definition 
a{f) is anti-linear in /), and they satisfy the canonical commutation relations 

Mf),a*(g)] = (f,g) L i m , [a(f),a(g)} = [a*(f),a*(g)} = 0. (4.6) 
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For every / E L (R ), we introduce the self adjoint operator 

<P(f)=a*(f) + a(f). 
We will also make use of operator valued distributions a* and a x (x E IR 3 ), denned so that 

a*(/)= [dxf(x)a x 



a (f) = / dx/(x)a x . 

for every / E L 2 (R 3 ). The canonical commutation relations take the form 

[a x , a*] = 8{x - y) [a x ,a y ] = [< , a*] = . 
The number of particle operator, expressed through the distributions a x ,a*, is given by 

N = l dx a*a x . 



The following lemma provides some useful bounds to control creation and annihilation operators 
in terms of the number of particle operator M. 

Lemma 4.2. Let f E L 2 (R 3 ). Then 

IK/MI < ll/H HA^H 

||a*(/MI< 11/11 II (AA + 1) 1/2 ^II 
||<A(/)V||<2||/||||(AA + 1) 1 / 2 ^||. 

Proof. The last inequality clearly follows from the first two. To prove the first bound we note that 

\\a(f)il>\\<Jdx\f(x)\\M\\ < (Jdx\f(x)\ 2 ^) (Jdx\\a x iP\A 

= 11/11 IIAT^II. 



The second estimate follows by the canonical commutation relations (|4.6p because 

IK Wll 2 = ftM/)a*(/ty> = (i>,a*(f)a(m + ||/|| 2 |H| 2 

= \wm\ 2 + ii/ii 2 m 2 < ii/ii 2 (w l/2 n + u\\ 2 ) = ii/ii 2 ii (M+l) l ' 2 n 2 



a 



Given ip E J 7 , we define the one-particle density 7^ associated with ip as the operator on L' 



,Wi 1 



with kernel given by 

By definition, j\ is a positive trace class operator on L 2 (R 3 ) with Tr j\' = 1. For every iV-particle 
state with wave function i/jn E L 2 (R 3Af ) (described on T by the sequence {0, 0, ... , ipN,0, 0, . . . }) it 
is simple to see that this definition is equivalent to the definition (jl.4p . 
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We define the Hamiltonian H N on T by {H N ^ n) = H^ } ^ n) , with 

n i n 

jf=l i<3 

Using the distributions a x , a*, TCn can be rewritten as 

H N = dxV x a* x X7 x a x + — dxdyV(x - y)a* x a y a y a x . (4.8) 

By definition the Hamiltonian T~Ln leaves sectors of J- with a fixed number of particles invariant. 
Moreover, it is clear that on the TV-particle sector, T~Lm agrees with the Hamiltonian Hn (the subscript 
N in TCn is a reminder of the scaling factor 1/N in front of the potential energy). We will study 
the dynamics generated by the operator TLn- In particular we will consider the time evolution of 
coherent states, which we introduce next. 

For / £ L 2 (]R 3 ), we define the Weyl-operator 

W(f) = exp (a* (/) - a(f)) = exp (J dx (f(x)a x - J(x)a x )^ . 
Then the coherent state ip(f) £ -F with one-particle wave function / is defined by 

Notice that 

m = w(m = e -" /l|2/2 E ^F n = e "" /l|2/2 E 4t f m > ( 4 - 9 ) 



n>0 n>0 

where Z®" - indicates the Fock-vector {0, . . . , 0, /® n , 0, . . . }. This follows from 

exp(o*(/) - o(/)) = e-ll/H 2 /2 exp(a*(/)) exp(-a(/)) 

which is a consequence of the fact that the commutator [a(f),a*(f)] = ||/|| 2 commutes with a(/) and 
a*(f). From (|4,9p . it follows that coherent states are superpositions of states with different number 
of particles (the probability of having n particles in ip(f) is given by e"' " \\f\\ 2n /n\). 

In the following lemma we collect some important and well known properties of Weyl operators 
and coherent states. 



Lemma 4.3. Let f,g € L 2 ( 

i) The Weyl operator satisfy the relations 

W(f)W(g) = W(g)W(f)e- 2iIm{f ' 9) = W(f + g) e - iIm{f > 9) . 

ii) W{f) is a unitary operator and 

w(fy = w(f)- i = w(-f). 

Hi) For every i£l 3 , we have 

W*(f)a x W(f) = a x + f(x), and W*(f)a* x W(f) = a* + J(x) 
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iv) From Hi) it follows that coherent states are eigenvectors of annihilation operators 

aMf) = f(x'Mf) => a(g)ip(f) = (g, f) L ^tf) . 

v) The expectation of the number of particles in the coherent state ip(f) is given by \\f\\ 2 , that is 

(iP(f),Mi>(f)) = 



Also the variance of the number of particles in ip(f) is given by \\f\\ 2 (the distribution of M is 
Poisson), that is 

(V(/),AT 2 ^(/)) - (V(/),A/W)) 2 = 



vi) Coherent states are normalized but not orthogonal to each other. In fact 

(i>(f),^(9)) = e-Hll/ll 2 +N! 2 -2(/, S )) ^ Mf),il>(g))\ = e-51!/-^" 2 . 

4.2 Time Evolution of Coherent States 

Next we study the dynamics of coherent states with expected number of particles N in the limit 
N — ► oo. We choose the initial data 



tp{VN<f) = W(VN<p)n for if G # X (R 3 ) with \\<p\\ = 1 
and we consider its time evolution ip(N,t) = e~ t ' l ~ iNt ip(\'iVip) with the Hamiltonian TC^ defined in 



Theorem 4.4. Suppose that there exists D > such that the operator inequality 

V\x) < D{1 - A x ) (4.10) 

holds true. Let T Nt be the one-particle marginal associated with if)(N, t) = e~ i ' HNt W \y/Ntp)Q (as 
defined in KJTfy). Then there exist constants C, K > (only depending on the H 1 -norm of if and on 
the constant D appearing in {J^.IO^ ) such that 

r$ t -\f t )(f t \ <%e K W (4.11) 



Tr 

for all t G R. 



N 



(4.12) 



We explain next the main steps in the proof of Theorem 14.41 By (|4.7[) , the kernel of TjJ t is given 
by 

r^foy) ~ ^,W*(VNf)e iHNt a* y a x e- iHNt W(^/Nf)^ 

= ft(x)f t (y) + ^ (O, W*(VNif)e m ^(a x - VNf t (x)) e - m ^W(VNif)n) 

+ ^0- (n,W*{VNif)e mNt {a* v - VNf t (y))e- inNt W{VNif)n) 
\/N \ ' y / 

+ j- (n,W*(VNf)e iHNt (a* y - VNf t (y))(a x - VNf t {x))e- mNt W{VNf)VL 
It was observed by Hepp in [21] (see also Eqs. (1.1 7)- (1.28) in [20j) that 

W*(VN(ps) e m ^- s \a x - VNif t (x))e~ mN ^-^W(VNif s ) =U N (t;sY a x U N (t;s) 

= U N (s;t) a x U N (t; s) 
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(4.13) 



where the unitary evolution Ujy(t; s) is determined by the equation 

id t U N {t;s) = C N (t)U N (t;s) and U N (s;s) = l (4.14) 

with the generator 

£jv(i) = / dxV x a* x V x a x + dx (V * \(ft\ 2 ) (x)a* x a x + / dxdyV(x - y)~^l(x)ipt{y)a y a x 

1 

2 



+ - I dxdy V(x - y) (ip t (x)ip t (y)a x a* y + ip t (x)<p t (y)a x a y ) 



+ -7= J dxdyV(x - y) a* ((pt(y)a>l + <Pt(y)a>y) 



(4.15) 



+ 2iV / dxdy V ^ X ~ ^ a * a *y a y a 

It follows from (|4TT2|) that 



r^ t (x,y) - <pt(x)lp t (y) = i <n,Z^(t;0)*^o x W2v(t;0)n) 



+ ( f^-(n,U N (t;0)*a* y U N (t;0)n) (4.I6) 

+ ^(ft,^(i;0)*a^ iv (i;0)tt> . 
v iv 



In order to produce another decaying factor 1/vJV in the last two term on the r.h.s. of the last 
equation, we compare the evolution Wjv(£;0) with another evolution ZYjv(£;0) defined through the 
equation 

id t U N {t]s) = C N (t)U N (t;s) with z7jv(s;s) = 1 (4.17) 



and 



£jv(£) = / dx V x a*V x a x + dx (V * \(p t \ 2 ) {x)a* x a x + / dxdy V(s - y)(p t (x)ip t (y)a* y a x 



+ - J dxdyV(x-y)((pt(x)(pt(y)a x a* y + pt(x)pt(y)a x a y ) 



+ 2N dxd y V ( x - y) a la* y a y a a 
From (|4.16p we find 

T { l\{x;y) - ip t {x)Tp t {y) 



^{n,U N (t;0)*a* y a x UN(t;0)n) 
"" " f, Sl t U N (t\G)*<% (u N (t;0) -U N (t;0)\ fi) + (0, (Z4r(t;0)* - Z^(t;0)*) a* y U N (t;0)n 



N 
+ ^M((0.// v (/;U)'V/. r (/Y.v(/:()!-//v(/:()))0) + ((l(/Y.v(/:()) i -U v </:<)) ) «. ,/Y v ( /: ())<> ) ) . 



(4.18) 
because 



tt,UN(t;0)*a y U N (t;0)n) = ( n,U N (t; 0)* a* U N (t; 0)n) = 0. 
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This follows from the observation that, although the evolution Wjv(£) does not preserve the number 
of particles, it preserves the parity (it commutes with (— 1)-™). From (|4.18p . it easily follows that 



,(i) 



1 



r$, - \<Pt)(<Pt\ < T7 (U N (t;0)n,MU N (t;0)Q) 



hs ~ N 



\\(U N (t; 0) - U N (t; 0))O|| || (AA + l) l l 2 U N (t- 0)O|| 
||(Wiv(t; 0) - Wjv(t; 0))O|| || (AA + l) 1/2 ^iv(t; 0)O|| 



(4.19) 



To bound the r.h.s. of (|4.19p . we need to compare the dynamics Wjv(i;0) and Uj^(t;0), and to 
control the growth of the number of particle N with respect to the fluctuation dynamics Z//jv(£;0) 
and Wjv(i; 0). We show, first of all, that 



(U N (t;0)n,NU N (t;0)n) <Ce x|t| . 
To prove this bound, we compute the time derivative 

— (U N (t; 0)O, (A/" + 1) U N (t; 0)O) = {U N (t;0)n, [C N (t),J^]U N (t;0)n) 
at 



(4.20) 



(^;0)O,[£,v(i),A^(t;0)O) 

21m / dxdyF(x - y) Vt {x) n {y){U N {t; 0)O, [o*oJ,JS/]Z<v(*; 0)O) 

41m / dxdyV(x - y)ipt(x)ip t {y)(UN(t; 0)O, a*a*^v(t; 0)O) . 



Thus, from Lemma 14.2} we obtain 

A(Z4r(i;0)O, (AT +l)U N (t; 0)O) 

dt 

< 4 fdx\cp t (x)\ \\aJJ N (jb;0)Sl\\ \\a*(V(x - .)cp t )U N (t; 0)n\\ 

< 4sup ||y(x - .)^|| ||(AA + l) 1/2 ^v(i; 0)O|| 2 

< C {U N (t; 0)O, (AA + l)U N (t; 0)O) , 
where we used the fact that 



(4.21) 



\\V(x-.)tp t \ 



dy|y(x-2/)| 2 |^(y)| 2 <C||^||^<C||^| 



H l 



because of the assumption (|4.10p . From (|4.2ip . we obtain (|4.20p applying Gronwall's Lemma. 

Making use of (|4.20p (and of an analogous bound for the growth of the expectation of M 4 w.r.t. 
the evolution U]y{t; 0); see [28] [Lemma 3.7]), we can derive the bound 



U N (t;0) -U N (t;0)) O 



< 



C 



,K\t\ 



N 



(4.22) 



for the difference between the two time evolutions Un{1\ 0) and ^at(£; 0) (note that, at least formally, 
the difference between the two generators Cjq{t) and £jv(£) is a term of the order N~ 1 ' 2 ; this explains 
the decay in iV on the r.h.s. of (|4.22p ). 
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In [2H [20] the time evolution U(t; s) was proven to converge, as N — > oo, to a limiting dynamics 
Moo(t] s) denned by 

id-tUooit; s) = Coo(t)Uoo(t', s ) and Uoo{s; s) = 1 

with generator 

Coo(t) = / dx V x a* x X7 x a x + / dx (V * |</2 t | 2 ) (x) a* x a x + / da?dyF(a? - y)Tpi(x)<p t {y)a y a x 

+ 2 / dxdy ^( x ~ y ) (^fa)^ (2/)°& a » + 1 Pi{x)Tpi(y)a x a y ) . 

In this sense, Hepp (in [21], for smooth potentials) and Ginibre-Velo (in [20] . for singular potentials) 
were able to identify the limiting time evolution of the fluctuations around the Hartree dynamics. 
Analogously to (|4.22p . the bound (|4.20p can also be used to show that, for a dense set of ^ E J 7 , 
there exists constants C,K such that \\(UN(t;s) — Uoo(t; s))^\\ < CN~ 1 l 2 e K \ t ~ s]{ , giving therefore a 
quantitative control on the convergence established in |21| |2"U] . 

Note, however, that the convergence oiUN{t;s) to the limiting dynamics Uoo(t;s) is still not 
enough to obtain estimates on the difference between r^- t and \(p t )(ipt\- In fact, to reach this goal, 
we still need, by (|4.19p . to control the growth of J\f with respect to the time evolution Z^at(£; s). We 
are going to prove that 

(U N (t;0)Q,AfU N (t;0)n) < Ce KW . (4.23) 

Inserting (J4T20J) . P~22j) and (J433|) on the r.h.s. of P~l~9|) . it follows immediately that 



,d) 



,e 



K\t\ 



r£J-kft><<ftl ns < c —> ( 4 - 24 ) 

which implies the claim (|4,lip . In fact, since \ipt){ipt\ is a rank one projection, the operator £7 = 
>Nt ~ l < /'*)(v 7 t| has at most one negative eigenvalue. Noticing that Tr<57 = 0, it follows that £7 has 
a negative eigenvalue, and that the negative eigenvalue must equal, in absolute value, the sum of all 
its positive eigenvalues. Therefore, the trace norm of £7 is twice as large as the operator norm of £7. 
Since the operator norm is always controlled by the Hilbert Schmidt norm, we obtain (|4.1ip (this 
nice argument was pointed out to us by R. Seiringer). 

The proof of (J4.23P is much more involved than the proof of the analogous bound (|4.20p . This 
is a consequence of the presence, in the generator C^(t), of terms which are cubic in the creation 
and annihilation operators (these terms are absent from Cft^t)). Because of these terms, also the 
commutator [£tv(£),AA] is cubic in creation and annihilation operators, and thus its expectation (in 
absolute value) cannot be controlled by the expectation of M . For this reason, to prove (|4.23p . we 
have to introduce yet another approximate dynamics Wn(P, s), defined by 

id t W N (t; s) = JC N (t)W N (t; s), with W N {s; s) = 1 

and with generator 

fCN{t) = / dxV x a* x V x a x + / dx (V * \<p t \ 2 ) {x) a* x a x + / dxdy V(x - y)Tpi{x)ip t {y)a* y a x 



+ 2 / dxd V v ( x ~ v) {<Pt(x)<p t (y)alal + (p t (x)(p t (y)a x a y ) 

+ 4= / dxd V V(x - y) a* x {<p t (y)xW < N)a* y + Tp- t (y)a yX (M < N)) 



(4.25) 



+ 2iV / dxdy V ( x ~ y " s a *x a *y a y a x ■ 
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Observe, that the generator K,jy(t) has exactly the same form as the generator £j\r(i); the only 
difference is the presence of a cutoff in the number of particles J\f inserted in the cubic term. Thanks 
to the cutoff in AT and to the factor N^ 1 ' 2 in front of the cubic term in ICj^(t), we can prove, making 
use of a Gronwall-type argument, that 

(W N (t;s)n,AfW N (t;s)n) <Ce K ^- s \. (4.26) 

Actually, it is simple to see that the last inequality can be improved to 

(W N (t; s) n,M j W N (t; s)0) < Cj e K ^- s] . (4.27) 

for every j£i and for appropriate j-dependent constants Cj and Kj. To obtain ()4.23p . we still have 
to compare the dynamics Z^7v(t; s) and Wat(£; s). To this end, we first show weak a-priori bounds of 
the form 

(U N (t;s)i;,Af j U N (t;s)i;} < C {t/>, {N + N + 1) J » (4.28) 

for every ip G T and for j = 1, 2, 3 (these bounds hold uniformly in t, s £ M and can be proven using 
the very definition of the unitary group U]y(t; s); see [28] [Lemma 3.6]). Using (|4.28p . we find 

(U N (t; 0)n,AT{U N (t; 0) - W N (t; 0)) 0) 

= (U N (t; 0)n,ATU N (t; 0) (1 - U N (t; 0)*W N (t; 0)) O) 

= -i f ds (U N {t;0)n,NU N {t;s) (£ N (s) - lC N (s))W N {s;0)n) 
Jo 

r-t 



i 

N ./„ 



ds / dxdyV(x — y) 



K (U N (t; 0)Q,MH N (t; s)a* x (^ t {y)a yX {N > N) + <pt(y)x(tf > N K) a x W N (s; 0)O) 
ds / dx (a x U N (t; s)*MU N (t; 0)0, a(V(x - .)<p t )x(M > N)a x W N (s; 0)0) 



NJ 

t 
ds I dx{a x U N {t- s)*MU N {t] 0)O, X {M > N)a*{V{x - .)ip t )a x W N (s; 0)0) . 



<N Jo 
Hence 

(U N (t; 0)n,Af(U N (t; 0) - W N (t; 0)) O) 

< -L [ ds [ dx\\a x U N (t;s)*NK N (t;0)n\\\\a(V(x - .)ipt)a xX (N > N +l)W N {s;0)n\\ 
y/N Jo J 



+ -^= [ ds [ dx\\a x U N (t;syMU N (t;0)n\\\\a*(V(x - .)^ t )a xX (M > N)W N (s;0)n\\ 
v N Jo J 

< 2su P*H^--)^H f <± s \\M l / 2 u N {t-syNu N (t- o)vil \\M X (M > N)w N (a;o)ii>\\ . 

vN Jo 

Therefore, using the inequality X (N > N) < (N/N) 2 and applying (J4.27P (with j = 4) and (|4.28p 
(first with j = 1, and then with j = 3) we can bound the two norms in the s-integral. It follows that 



(U N {t- 0)n,M(U N (t; 0) - W N (t; 0)) 0) 
which, combined with (|4.26p . implies (|4.23l) . 
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<Ce 



Kt 



4.3 Time-evolution of Factorized Initial Data 

To prove Theorem 14.11 we express the factorized initial data as a linear combination of coherent 
states. Using the properties listed in Lemma 1431 it is simple to check that 

{o, o, . . . , o, <p® N , o, o, . . . } = {a * { ^} N n = d N r — e im w(e- ie VNv)n 

VNl Jo 27T 



with the constant 



d - = w^m * Nl/A 



The kernel of the one-particle density jj^ t associated with the solution of the Schrodinger equation 
{0, . . . , 0, e~ iHtft (p e > N ,0, . . . } is thus given by (see g2D) 

= f- f' '^ t '^e-' e '"e«'^(W(e-^VN V )n,a' v (t)a x (t)W{e-«''VN^n) 



where we introduced the notation a x (t) = e l ^ Nt a x e l ^- Nt . Next, we expand 

72 [ 2 * d9i f 2lT 
L ^ Jo 



x (a x (t)-e- i02 VN<f t (x)] W (e~ i92 VN ^)ft\ 



^N Jo 2tt J 2tt 



x (w(e~ i0l VNip)n, (a x (t) - e- ie2 VNip t (x)) W(e~ ie2 VNip)n\ 



y/N Jo 2vr J 2vr 



W(e- i01 VNip)n, (a*(t) - e i6l VNlp t (y)\ W (e^ 2 VN (p)^ 



^de 1 f^de 2 m 



W{e~ iei VN<p)n, W(e- i62 VN(f)n\ . 



(4.29) 



Since 



d N 



F ^.e i0 ( N -^W(e- i6 ^/N<p)n = d N e~ N ' 2 f] ( F ^ e "<"-wA N^^^-Sl 



3=0 

e -N/2 N (N-l)/2 (a * (¥j)) iV-l 
«iV , — —. 1 1 

y/N-V. N - 1! 
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we find that 

_ = d N f * dOi f w d02 e - iei N e ie2N 
' N J 2vr J Q 2ir 



x (W(e-™WN<p)n, (a*(t) 



N<p t (y)\ (a x (t) 



-102 



Nip t {x))W(e- ty2 VN^)Q 



f ^M [ n ^1 e^ N (pW-V, (a x (t) - e-^VNip t (x)) W( e -^VN<p)n 



+ 



d N (f t (x) 



N Jo 2vr 

2vr 



N 



I ' I 1 e ~ l9lN (wie-^VN^a, (a;(t) - e*VmW) V*"^ 



(4.30) 



and thus 



(i) 



-y { N > t (x;y)-(pt(x)p t (y) 



2n d9o 



j2 f-2ir j n 

-at/o ^J ^hy^N^m\\\\a x K 2 (f,o)n\\ 
*,M*)| /- 2 '^ ||arf(t . 0)si|| 



A^ 



where U^(t;s) is defined as (|4.14p . with ip t replaced by e ipt (we are using, here, the fact that, 
although the Hartree equation is nonlinear, e ipt is always a solution if ipt is). Since cLn — N 1 ' 4 , it 
follows from the bound (|4.23|) for the growth of the expectation of N with respect to the fluctuation 
evolution U%(t; s) (the bound clearly holds uniformly in 6) that 



(!) I W I 



< 



c 



Jit 



hs ~ N 1 / 4 
and therefore (using the argument presented after (|4.24|) that 



Tr 



(!) i \/ 



< 



C 

A7I74' 



„i« 



(4.31) 



To improve the decay in A" on the r.h.s. of (|4.3ip from A^ -1 ' 4 to N~ 1 ' 2 (as claimed in The- 
orem 14. ip . it is necessary to study the second and third error terms on the r.h.s. of (|4.30p more 
precisely; for the details, see [28] [Lemma 4.2]. 



A Non-Standard Sobolev- and Poincare Inequalities 

In this section, we collect some non-standard Sobolev- and Poincare-type inequalities which are very 
useful when dealing with singular potentials. 

Lemma A.l (Sobolev-type inequalities). Let -0 £ L 2 (M 6 , dzida^)- IfVE. L 3 / 2 (M 3 ), we have 

(V, V(xi - * 2 )V> < C||y|| 3/2 (V, (1 - AOV) . (A.l) 

IfV eL^R 3 ), then 



(V>, V(xi - x 2 )^) < C\\V\\i (V, (1 - Ai)(l - A 2 )i>) 



(A.2) 
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The first bound follows from a Holder inequality followed by a standard Sobolev inequality (in 
the variable xi, with fixed x 2 ). The proof of (|A.2p can be obtained through the same arguments 
used in the proof of the next Poincare-type inequality (see [T5]). 

Lemma A. 2 (Poincare-type inequality). Suppose that h E L 1 (M 3 ) is a probability density with 
J dx \x\ 1 ' 2 h(x) < oo. For a > 0, let h a (x) = a~ 3 h(x/a). Then we have, for every < k < 1/2, 

\(<p, (ha(x! - x 2 ) - 8( Xl - x 2 )) ^)| < Ca K (<p, (1 - Ai)(l - & 2 )<p) X/2 (il>, (1 - Ai)(l - A 2 )^) 1/2 . 
Proof. We rewrite the inner product in Fourier space. 
(<p, [h a (xi - x 2 ) - S(xi - x 2 )Ji(j) 

dpidp 2 dqidq 2 dx6(pi+p 2 -q 1 -q 2 )(p(pi,p 2 ) : ip(qi,q 2 )h(x) ^ e m(pi_(?l) '' - i 
Using that \ e ^{vi-qi)-x _ i| < a K |x| K |pi - gi| K , we obtain 
(ip, (h a (xi - x 2 )-5(xi - x 2 ))i(j) 

x / dpidp 2 dq 1 dq 2 5(p 1 + p 2 - qi - q 2 ) (\pi\ K + \qi\ K ) \<p(pi,P2)\ |^(9l,92)| • 
We show how to control the term proportional to \pi\ K ; the other term can be handled similarly. 

(ip, (h a (xi - x 2 ) - 5(x 1 - x 2 )) tp) 

/l IK/1 

dpidp 2 dqidq 2 8 (pi + p 2 - qi - q 2 ) 1 



bi r(l+pj) (1 -" )/2 (l+pj) 1/2 i^ v, 

(l + ^)V2(l + g |)i/ 2 l^'^l 



< ^M^±^L^qi,q 2) \ 



(l+p2)(l-*)/2( 1+p 2)l/2 

(i+^)(i + p!) x1/2 

+ 9 2 2 ) 



< Ca K ( [ dpidp 2 dqidq 2 8(pi +p 2 -q x - g 2 ) ( * t^j T^ lfi(Pi,P2)| 2 

V (! + 9i)(l + 92) 

H H H H A^ x , (l + g 2 )(l + q 2 ) ? ^ 1/2 

dpidpadgidga «5(Pi +P2 - 9i - ft) (1 + p2) i- K(1 + p 2) hKgi.gs)l 

< Ca 1/2 <^, (1 - Ai)(l - A 2 ) V 9) 1 / 2 (V, (1 - Ai)(l - A 2 )V) 1/2 

x (?/ d9 d +? 2 )(i 1 + ( P - ? ) 2 )) ( s r/ dp < 



l+p 2 )(l + (g-p) 2 ) 



2M-K 



1/2 



The claim follows because 



sup / dp 



(l+pi)(l + ( q -pf) 



2U-K 



<C 



(A.3) 



for all k < 1/2. To prove (|A.3jl we consider the three regions \p\ > 2\q\, \q\/2 < \p\ < 2\q\ and 
\p\ < \q\/2 separately. Since \p — q\ > \p\/2 for \p\ > 2\q\, it follows that 



dp 



< 



P |>2| g |(l+P 2 )(l + ('Z-P) 2 ) 1 - K ~J\ P \>2\ q \ ( 1 + pL 



dp 



1—K 



<C 



dp 



(l+P 2 ) 



2 \2-K 



< OO 
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for k < 1/2, uniformly in q. For \p\ < \q\/2, we use the fact that \q — p\ > \q\/2, and we obtain 

f dp C f dp C\q\ 

4|<M/2 (1 + P 2 )(l + (q ~ PW- K ~ (1 + Q 2 V- K 4|<|,|/ 2 1 + P 2 - (1 + q 2 ) 1 - 

which is bounded uniformly in q. Finally, in the region \q\/2 < \p\ < 2\q\, we use that 



iM/2<| P |<2lal(l+P 2 )(l + (9-p) 2 ) 1 - K " a+q 2 )J\ P \ 



dp , J9| 2k+1 ^ 

< C TT < OO 



m/ 2 <h<2m(i+p 2 )(i + (9-p) 2 ) 1 - k " (i + ? 2 )7h<3m(i+p 2 ) 1 - k - 1 + 9 2 

uniformly in g 6 I 3 , for all n < 1/2. D 



In the approach developed in [15] for the case of large interaction potential we can only prove 
weaker estimates on the solution ipu,t of the Schrodinger equation. As discussed in Section 13.31 we 
can only prove that 

(W* N , (i>j) i>N,t, ((V, • V,) 2 - A, - Aj + 1) w* NM ^ N , t ) < C 

uniformly in N and t. For this reason, we need estimates which only require the boundedness of the 
expectation of this particular combination of derivatives. The next lemma gives a Sobolev inequality 
of this type. 

Lemma A. 3. Suppose V G L 1 (IR 3 ). Then 

\(<p, V(xi - x 2 )il))\ < C\\V\\x (V, ((Vi • V 2 ) 2 - Ai - A 2 + 1) V) 1/2 

x (ip, ((Vi • V 2 ) 2 - A! - A 2 + 1) vj) 1 / 2 

for every tp, (p G L 2 (R 6 ,dxidx 2 ). 

Proof. Switching to Fourier space, we find 

(tp,V(xi ~x 2 )ip) = / dpidp 2 dqidq 2 <f(pi,P2)4>(qi,q2)V(qi -p\)8(p\ + p 2 - qi -92)- 

Therefore, by a weighted Schwarz inequality, 

{tp,V(xi -x 2 )ip) 

< F 00 / dpidp 2 d^idg 2 -^ — 3— — — — - MPi,p 2 )\ d(pi +P2 - ?l - <?2j 

/ dpidp 2 dgidg 2 -^ — 2 \ \ \^{qi,q 2 )\ 2 5{pi + p 2 - qi - q 2 ) ) 

J {pi-P2r+Pi+p 2 + i J 



~ mh { S 7l dq (q-(p-q)) 2 + l 2 + (p-q) 2 + l 

x (i/>, ((Vi • V 2 ) 2 - Ai - A 2 + 1) V> V2 (<p, ((Vi • V 2 ) 2 - Ax - A 2 + 1) vf 2 



The lemma follows from 



SUP ' dQ (q-{p-q)) 2 +l 2 + {p-q) 2 + l < °° " ^ 



53 



To prove (|A.4|) . we write 

d Q T-^Z _nno | o | 7- ^^^ = dq 



(q-(p-q)) 2 + q 2 + (p-q) 2 + l~~J\ q - ¥>lp \ n q _Ef_EL\\ q 2 + {p _ q? + l 



+ / dg 2 



(A.5) 



The first term on the r.h.s. of the last equation is bounded by 

M ~, : m> ^ / d <? 



2 „2\ 2 . o . / no . , ./l„_£l^l„l AL_2| 4 



16 /" , 

< "5" / d? — d ■ - < oo 



k-fl>H ((g_|) 2 _l£) +(Z 2 + (p _ (z) 2 + 1 7|«-f|>[p| ^|9-f| +1 

1 

uniformly in p € R 3 . As for the second term on the r.h.s. of (1A,5|) . we observe that 

L- iq - - -»' 



«-5Kbl ((,_|) 2 _£j- + ,2 + (/) _ g) 2 + | 

dx - 



/•bl r 2 

4-7T / dr g 

■A) (r 2 -M!) +2r 2 + ^ + l 

H/2 j 



< C\p\ 2 / dr 

J— Ie.I/2 



Pi/2 r 2 (r + | p |)2 + ^ + M) 2 + M! + l 

'bl/2 1 /• 1 
<C I dr — 7 <C/dr — < oo, 

bl/2 



r 2 + 1 ,/k r 2 + 1 



uniformly in p. □ 
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